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In  this  project  we  studied  topological  phases  at  transition  metal  oxide  interfaces.  Over  the  course  of  the  5-year  project,  a  fairly  detailed 
understanding  of  these  systems  was  acquired.  By  using  a  combination  of  model  Hamiltonian  studies  and  first-principles  calculations,  we 
made  a  number  of  predictions  for  specific  material  systems  that  might  exhibit  topological  phases,  with  the  most  promising  candidate  being 
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ARO  PECASE  Research  Summary 
Gregory  A.  Fiete 
University  of  Texas  at  Austin 


The  title  of  my  PECASE  project  was  "Novel  quantum  phases  at  interfaces”.  The  main  focus 
of  this  project  was  the  exploration  of  novel  quantum  phenomena  in  correlated  materials  at 
the  interface  of  transition  metal  oxides,  hut  over  the  5-year  duration  of  the  award  many 
closely  related  projects  were  undertaken  that  could  help  to  inform  the  understanding  of 
quantum  behavior  at  interfaces.  In  particular,  the  study  of  topological  phases  featured 
prominently  throughout;  1  believe  the  work  funded  here  has  had  a  broader  impact  on  the 
field  of  condensed  matter  physics  by  bringing  others  into  this  line  of  inquiry  and  by  laying 
some  of  the  initial  (and  current)  expectations  for  how  likely  one  is  to  find  certain  correlated 
electronic  phases  in  the  interfaces  and  heterostructures  we  studied.  Eor  the  DoD,  and  Army 
in  particular,  1  hope  the  research  summarized  here  will  ultimately  be  useful  for  next 
generation  technologies,  particularly  in  the  areas  of  communication  and  sensing. 

Approximately  50  papers — and  a  similar  number  of  invited  talks — resulted  from  this  award 
(including  works  in  Nature  Physics,  Nature  Materials,  Science,  and  Physical  Review  Letters), 
so  it  is  not  possible  to  summarize  all  of  the  results  in  detail.  The  most  important  results 
were  conveyed  in  real-time  to  the  program  manager  for  this  project.  Dr.  Marc  Ulrich.  In 
order  to  best  facilitate  the  internalization  of  some  of  the  most  important  results,  1  have 
provided  a  "lay"  summary  in  the  text  that  immediately  follows,  and  1  have  also  attached  3 
review  articles  (either  in  draft  form,  or  the  final  published  form,  if  available).  Together,  they 
provide  an  excellent  summary,  with  references,  of  the  most  important  technical 
accomplishments  of  this  project.  The  first,  "Topological  Phases  in  Oxide  Heterostructures 
with  Light  and  Heavy  Transition  metal  ions"  highlights  work  at  the  core  of  this  project.  The 
second,  "Eractionalized  Topological  Phases",  highlights  some  of  the  new  ground  broken  on 
the  fundamental  influence  of  interactions  on  topological  phases.  The  third,  "Topological 
Insulators  and  Quantum  Spin  Liquids",  highlights  deep  connections  between  frustrated 
magnetic  systems  with  no  magnetic  order  and  topological  insulators.  In  addition  to  the 
technical  accomplishments,  many  students  and  postdocs  were  trained,  and  countless 
informal  discussions  were  held  with  professional  colleagues  in  the  dissemination  of  these 
results. 


Some  Specific  Accomplishments  from  the  Last  Few  Years 

Recently,  1  have  focused  my  group's  research  efforts  on  two  main  topics  at  the  forefront  of 
research  in  condensed  matter  physics: 

(i)  Topological  Insulators 

(ii)  Quantum  Spin  Liquids 

These  names  may  not  be  especially  meaningful  to  scientists  outside  of  condensed  matter 
physics,  so  1  will  briefly  summarize  what  they  are,  as  well  as  a  few  key  contributions  of  my 
group  to  each  of  these  areas. 

Topological  Insulators 

Most  materials  fall  into  one  of  two  categories — a  metal  (conductor)  or  an  insulator. 
(Semiconductors  are  not  "half-conductors",  but  rather  insulators  that  are  in  a  certain  sense 
"weakly"  insulating.)  For  about  half  a  century,  physicists  thought  they  understood 
insulators  of  the  semiconducting  variety  rather  well.  We  have  designed  the  transistor 


around  their  properties,  and  most  electronic  devices  incorporate  them  in  some  way.  Such 
insulators  are  commonplace. 


It  thus  came  as  a  major  surprise  in  the  first  decade  of  the  21^^  century  that  physicists  had 
missed  something  huge  in  this  class  of  materials:  The  full  breadth  of  relativistic  effects.  The 
theory  of  Special  Relativity,  for  which  Einstein  is  famous,  plays  a  role  in  solids  when  the 
nuclei  are  heavy,  and  this  modifies  the  electronic  motion  in  highly  non-trivial  ways.  In  the 
language  of  solid  state  physics,  this  effect  is  called  "spin-orhit  coupling"  because  the  spin,  or 
magnetic  moment  of  the  electron,  is  coupled  to  its  motion  through  a  relativistic 
transformation  of  the  nuclear  electric  field  to  a  magnetic  field  in  the  reference  frame  of  the 
electron.  As  the  electron  moves  through  the  solid,  the  direction  of  its  spin  changes  in 
response  to  the  varying  magnetic  field  in  its  frame.  While  this  effect  has  been  known  for 
some  time,  it  was  previously  thought  have  a  very  weak  on  the  overall  physics  of  materials. 
Physicists  have  since  learned  that  spin-orbit  coupling  can  drive  a  phase  transition  from  an 
ordinary  insulator  to  a  "topological  insulator". 

A  topological  insulator  (of  the  most  common  variety — "Z2")  has  a  precise  mathematical 
description,  but  its  salient  property  is  that  it  is  insulating  in  the  bulk  and  conducting  on  its 
boundary  (surface).  Moreover,  its  conducting  surface  is  unlike  any  previously  known 
conductor,  and  it  is  predicted  to  have  a  number  of  potentially  novel  uses  in  electronic 
devices,  including  nearly  dissipationless  charge  and  spin  transport  at  room  temperature. 

In  a  rare  chapter  of  physics  history,  the  existence  of  the  topological  insulator  was  predicted 
theoretically  before  it  was  discovered  experimentally  (only  a  year  after  its  prediction).  The 
usual  order  in  condensed  matter  physics  is  experimental  discovery  first,  theoretical 
explanation  (if  any — high  temperature  superconductivity  still  does  not  have  one)  second. 
Today,  the  study  of  topological  insulators  is  still  dominated  by  theoretical  physicists,  and 
many  open  questions  remain,  particularly  involving  what  happens  in  materials  with  strong 
electron  correlations  and  additional  symmetries  (crystalline,  for  example)  that  may  enrich 
or  further  protect  topological  properties. 

The  issue  of  electron  interactions  lies  at  the  heart  of  my  group’s  effort  in  this  area,  and  more 
generally  in  all  the  problems  we  investigate  in  my  group.  Our  work  has  so  far  focused 
primarily  on  model  Hamiltonian  studies  in  which  we  predict  new  phases  of  matter  driven 
by  the  presence  of  both  strong  interactions  and  spin-orbit  coupling.  We  have  we  predicted 
a  new  phase  known  as  a  "weak  topological  Mott  insulator",  a  new  phase  known  as  "QSH*",  a 
new  phase  known  as  "TCMl”,  a  new  phase  known  as  "Tl*".  Each  of  these  phases  possesses 
some  unusual  properties  that  are  related  to  topological  properties  of  the  many-body  wave 
function  that  originate  in  strong  interactions  and  strong  spin-orbit  coupling.  These  phases 
cannot  exist  without  the  presence  of  significant  electron-electron  interactions,  and  in  this 
sense  are  truly  distinct  from  the  class  of  topological  phases  reported  in  bulk  Bi2Se3  and 
HgTe/CdTe  quantum  wells. 

Electronic  wave  functions  with  topological  properties  are  heavily  sought  after  because  they 
endow  usually  "fragile"  quantum  properties  with  orders  of  magnitude  improvements  in 
their  robustness.  Materials  with  topological  electronic  properties  can  open  the  door  to 
transformative  technologies  in  sensing  and  computing,  for  example. 

A  new  direction  undertaken  by  my  group  recently  is  the  theoretical  search  for  materials 
that  will  exhibit  topological  properties  predicted  in  model  Hamiltonian  studies  that 


incorporate  interactions  in  a  fundamental  way.  This  new  effort  will  bridge  theoretical 
studies  my  group  has  carried  out  with  the  latest  developments  in  materials  growth  [for 
example,  molecular  beam  epitaxy — used  to  grow  new  materials  atomic  layer-by-atomic 
layer).  Our  work  now  combines  state-of-the-art  numerical  methods,  such  as  DFT+DMFT 
and  DFT+Gutzwiller,  with  input  from  leading  experimental  groups  to  advance  technology  in 
this  direction.  In  addition,  we  are  also  studying  non-equilibrium  physics  in  models  [many 
relevant  to  topological  phases)  that  are  aimed  at  bringing  a  higher  degree  of  experimental 
relevance  than  many  of  those  studied  to  date,  including  the  "Floquet"  states. 

Quantum  Spin  Liquids 

Quantum  spin  liquids  are  a  rare  class  of  systems  with  local  magnetic  moments  or  "spins” 
which  are  not  macroscopically  magnetic  in  any  way — even  at  the  absolute  zero  of 
temperature.  In  effect,  the  magnetic  degrees  of  freedom  are  disordered  and  random,  much 
like  the  atomic  positions  in  a  liquid,  hence  the  name.  In  a  ferromagnetic  material  like  iron, 
most  of  the  magnetic  degrees  of  freedom  [spins)  orient  themselves  in  the  same  direction  to 
form  a  material  with  a  net  magnetization  at  room  temperature.  We  have  all  seen  how  this 
magnetization  can  be  used  to  pick  up  or  move  small  iron  filings,  for  example.  It  turns  out 
that  in  materials  with  local  spins,  the  spins  nearly  always  order  in  some  way — ferromagnetic, 
antiferromagnetic,  spiral,  etc.  at  some  finite  temperature.  This  is  because  of  the  presence  of 
a  characteristic  interaction  [exchange)  energy  scale  between  spins  [think  dipole-dipole 
interactions  from  electricity  and  magnetism,  for  simplicity)  that  will  become  important 
once  the  temperature  drops  below  this  energy  scale. 

However,  in  some  interacting  spin  systems  competing  energies  make  the  lowest  energy 
configuration  of  spins  difficult  to  "find".  This  feature  goes  under  the  general  name 
"frustration".  For  example,  consider  a  triangle  with  spins  on  each  vertex  [for  a  total  of 
three).  Suppose  the  energy  of  two  spins  is  minimized  if  they  are  oriented  in  exactly  the 
opposite  direction  [as  would  occur  in  an  ideal  antiferromagnet).  By  drawing  a  simple 
picture,  it  is  clear  that  it  is  not  possible  to  have  each  of  the  three  spins  oriented  opposite  to 
the  two  other  spins.  This  is  an  example  of  "frustration".  In  physical  systems  with 
frustration,  it  is  possible  that  the  system  is  never  able  to  find  a  configuration  that  minimizes 
the  total  energy  as  the  temperature  is  lowered,  and  it  fluctuates  [quantum  mechanically) 
between  many  different  configurations.  When  this  happens  at  zero  temperature,  the 
system  is  referred  to  as  a  quantum  spin  liquid. 

Why  would  one  be  interested  in  studying  quantum  spin  liquids?  The  study  of  quantum  spin 
liquids  is  vibrant  because  they  are  among  the  richest  [in  terms  of  phenomenology)  known 
interacting,  quantum  many-particle  systems.  As  a  result,  the  study  of  quantum  spin  liquids 
has  revealed  many  fundamental  aspects  of  quantum  mechanics.  It  is  probably  fair  to  say 
that  over  the  last  two  decades  a  fair  fraction  of  the  theoretical  advances  in  our  knowledge  of 
many-body  quantum  mechanics  is  related  to  quantum  spin  liquids  in  some  way.  There  are 
deep  connections  to  string  theory  via  the  technical  tool  of  conformal  field  theory, 
connections  to  high-energy  particle  physics  via  lattice  gauge  theories,  connections  to  the 
fractional  quantum  Hall  effect  via  fractionalized  quantum  numbers,  connections  to  quantum 
information  through  entanglement  studies,  and  promising  potential  for  quantum  computing 
via  non-Abelian  excitations  [notably  Majorana  fermions)  in  some  quantum  spin  liquids. 
Moreover,  quantum  spin  liquids  display  a  cornucopia  of  topological  properties. 

My  group  has  focused  on  exploring  the  topological  properties  of  the  many-body  wave 
functions  of  quantum  spin  liquids.  We  have  established  some  rather  general  statements 


and  powerful  relations  between  quantum  spin  liquids  and  topological  insulators,  unifying 
these  two  areas  of  research,  which  are  most  often  thought  of  as  disjoint,  it  is  one  concrete 
example  of  the  benefits  of  working  outside  the  traditional  boundaries  of  sub-fields  within 
condensed  matter  physics.  We  have  also  established  some  important  results  in  quantum 
spin  liquids,  such  as  that  disorder  in  the  Hamiltonian  can  help  stabilize  certain  types  of 
topological  phases,and  we  have  established  a  first  example  of  a  quantum  spin  liquid  phase. 


Summary  and  Future  Directions 

in  the  short  description  above  of  my  research,  i  left  out  a  number  of  other  areas  in  which  my 
group  has  made  contributions,  in  the  interest  of  brevity,  and  to  highlight  the  diversity  of 
our  research  efforts,  i  will  just  flag  the  work  on  the  entanglement  spectrum  of  strongly 
correlated  one-dimensional  systems  in  references. 

in  addition  to  the  basic  research  activity  described  here,  i  have  also  written  invited 
"Perspective”  articles  in  Science,  "News  &  Views"  in  Nature  Materials  and  invited  "Viewpoint" 
articles  in  Physics  intended  to  bring  important  developments  to  a  wide  audience,  i  have 
also  received  an  invitation  from  the  US  National  Academy  of  Sciences  to  address  the  Board  of 
Physics  and  Astronomy  in  a  45  minute  talk  on  the  potential  applications  of  topological 
insulators  in  devices:  http://sites.nationalacademies.org/BPA/BPA  067816.  Former  group 
members  have  gone  to  Harvard,  Stanford,  UC  Berkeley,  and  UIUC  among  other  places. 
Among  former  undergraduates,  three  have  received  NSF  Graduate  Research  Fellowships, 
and  one  was  an  APS  Apker  Award  Finalist. 

in  the  immediate  future,  my  group  will  be  devoting  a  significant  portion  of  its  research 
effort  to  studying  the  physics  of  transition  metal  oxides  with  heavy  transition  metals  that 
incorporate  both  strong  electron-electron  interactions  and  strong  spin-orbit  coupling — 
widely  acknowledged  as  key  ingredients  for  novel  quantum  phenomena  in  materials.  The 
transition  metal  oxides  have  already  given  us  high-temperature  superconductivity,  colossal 
magnetoresistance,  metal-insulator  transitions,  and  (through  molecular  beam  epitaxial 
growth)  can  be  used  to  produce  "materials  by  design"  with  layer-by-layer  growth,  it  is 
widely  thought  that  if  the  fundamental  physics  of  this  class  of  materials  can  be  understood 
and  controlled,  it  will  lead  to  a  new  generation  of  electronics  with  greater  functionality, 
sensitivity,  and  capability  than  currently  available. 

Finally,  there  is  no  need  to  restrict  research  effort  to  equilibrium  phenomena  in  this  or 
other  classes  of  materials.  The  rapid  development  of  laser  technology  in  the  last  decade 
combined  with  the  productive  cross-fertilization  between  cold  atom  physics  (where  non¬ 
equilibrium,  time-dependent  phenomena  are  readily  studied)  and  solid  state  physics  makes 
the  time  ripe  for  a  push  by  the  condensed  matter  community  to  address  non-equilibrium 
phenomena  in  correlated  materials.  My  current  funding  from  the  ARO  is  in  this  direction. 

We  have  already  developed  a  non-equilibrium  DMFT  code  that  presently  few  groups  in  the 
world  possess.  We  hope  that  we  will  be  able  to  settle  a  range  of  questions  of  both 
fundamental  interest  and  practical  relevance  to  real  experimental  systems. 


Topological  Phases  in  Oxide  Heterostructures  with  Light  and  Heavy 
Transition  Metal  Ions 

Gregory  A.  Fiete^'®’)  and  Andreas  Ruegg^ 

Department  of  Physics,  The  University  of  Texas  at  Austin,  Austin,  Texas  78712, 
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Theoretische  Physik,  ETH  Zurich,  CH-8093  Zurich,  Switzerland 
(Dated:  4  November  2014) 

Using  a  combination  of  density  functional  theory,  tight-binding  models,  and  Hartree-Fock  theory,  we  predict 
topological  phases  with  and  without  time-reversal  symmetry  breaking  in  oxide  heterostructures.  We  consider 
both  heterostructures  containing  light  transition  metal  ions,  and  those  containing  heavy  transition  metal 
ions.  We  find  the  (111)  growth  direction  naturally  leads  to  favorable  conditions  for  topological  phases  in  both 
perovskite  structures  and  pyrochlore  structures.  For  the  case  of  light  transition  metal  elements,  Hartree-Fock 
theory  predicts  the  spin-orbit  coupling  is  effectively  enhanced  by  on-site  multiple-orbital  interactions  and 
may  drive  the  system  through  a  topological  phase  transition,  while  heavy  elements  with  intrinsically  large 
spin-orbit  coupling  require  much  weaker,  or  even  vanishing  electron  interactions  to  bring  about  a  topological 
phase. 


I.  INTRODUCTION 

Topological  insulators  (TI)  have  now  been  studied  for 
nearly  a  decade  and  there  are  a  number  of  excellent  re¬ 
views  available,  from  both  an  experimental  and  theoret¬ 
ical  perspective. Aside  from  the  important  example 
of  SmBg,^^®  all  other  experimental  examples  of  topolog¬ 
ical  insulators  are  weakly  correlated.  The  most  impor¬ 
tant  reason  for  the  weak  correlation  effects  is  that  the 
bands  of  topological  insulators  near  the  Fermi  energy  are 
derived  from  s  and  p-type  orbitals,  which  are  rather  ex¬ 
tended.  By  contrast,  the  bands  about  the  Fermi  energy  in 
SmBg  are  predominantly  derived  from  /-orbitals,  which 
are  more  localized,  and  therefore  lead  to  flatter  bands 
and  enhanced  interaction  effects.® 

One  of  the  persistent  challenges  in  experimental  stud¬ 
ies  of  TI  is  the  problem  of  high  bulk  conductivity.'^  While 
some  progress  has  been  made  on  this  front  over  the  past 
few  years,  it  has  largely  been  incremental,  and  for  the 
most  part  has  been  focused  on  bismuth-based  TI.  An 
alternative  route  is  to  look  for  new  classes  of  strongly 
insulating  materials  that  might  support  topological  in¬ 
sulator  phases.  If  one  is  also  interested  in  studying  in¬ 
teraction  effects  that  could  possibly  take  one  beyond  the 
“band  theory”  description  of  topological  insulators,  ma¬ 
terials  with  orbitals  more  localized  than  the  s  and  p- 
orbitals  are  highly  desirable.  Transition  metal  oxides, 
which  typically  have  bands  derived  from  d-orbitals  close 
to  the  Fermi  energy,  are  excellent  candidates:  There  are 
a  large  number  of  insulating  oxides,  and  interaction  ef¬ 
fects  are  known  to  be  important  in  many  of  them-the 
high  temperature  cuprate  superconductors  serving  as  an 
excellent  example. Indeed,  there  have  been  a  number 
of  theoretical  proposals  for  strongly  correlated  topolog¬ 
ical  insulators  in  transition  metal  oxide  systems. A 
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recent  review  of  the  prospects  for  such  exotic  phases  in 
the  context  of  three  dimensional  iridium  (and  other  heavy 
transition  metals)  oxides  is  given  in  Ref.  [21]. 

In  this  article,  we  focus  on  topological  phases  in  thin 
film  (two-dimensional)  oxide  heterostructures  that  can 
be  described  within  the  band  theory  picture.  There  are 
three  prime  candidates:  (1)  The  time-reversal  invariant 
topological  insulator  characterized  by  a  single  Z2  invari¬ 
ant,  (2)  The  time-reversal  symmetry  broken  Chern  in¬ 
sulator  characterized  by  a  quantized  Chern  number  and 
quantized  Hall  conductance,  and  (3)  The  mirror  symme¬ 
try  protected  (with  respect  to  the  center  of  the  plane  of 
the  thin  film^^)  topological  crystalline  insulator  charac¬ 
terized  by  a  mirror  Chern  number. Combinations  of 
these  are  possible  as  well  for  a  “doubly  topological”  sys¬ 
tem,  though  clearly  (1)  and  (2)  are  mutually  exclusive. 
The  Chern  insulator  differs  from  the  quantum  Hall  in¬ 
sulator  in  that  the  former  has  time-reversal  symmetry 
spontaneously  broken  by  interactions  (that  is,  a  sponta¬ 
neous  magnetization  of  some  sort),  while  the  latter  has 
time-reversal  symmetry  broken  by  the  application  of  an 
external  magnetic  field.  As  a  result,  a  Chern  insulator 
requires  interactions  (because  there  is  no  spontaneous 
magnetism  without  electron-electron  interactions),  while 
a  (integer)  quantum  Hall  system  does  not.  A  topological 
crystalline  insulator  may  possess  time-reversal  symmetry 
or  have  it  broken  by  magnetic  order;  the  only  restriction 
is  that  any  magnetic  order  present  must  respect  the  mir¬ 
ror  symmetry. 

In  the  remainder  of  this  article,  we  will  consider  two 
transition  metal  oxide  structures,  the  perovskite  with  for¬ 
mula  ABO3  and  the  pyrochlore  oxide  A2B2O7,  where  A  is 
typically  a  rare-earth  element,  B  is  a  transition  metal  el¬ 
ement,  and  O  is  oxygen.  Two  specific  examples  we  study 
are  LaNiOa  and  ¥2^207,  which  are  both  materials  that 
have  been  grown  and  are  well  characterized  in  bulk  form. 
Our  new  angle  is  to  study  the  properties  of  thin  films  of 
these  materials  that  are  grown  along  the  (111)  crystalline 
axis.  We  find  the  physical  properties  of  these  films  are 
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rather  different  from  the  bulk.  Moreover,  there  does  not 
appear  to  be  an  obvious  way  to  infer  the  film  properties 
from  the  bulk.  As  a  result,  these  systems  appear  to  be 
excellent  candidates  for  exploring  novel  phenomena,  such 
as  topological  phases,  even  when  the  bulk  (non-thin  film) 
materials  possess  strikingly  different  properties,  such  as 
a  conducting  behavior.  We  are  thus  presented  with  the 
exciting  possibility  of  finding  “new  physics”  in  “old  ma¬ 
terials”  . 


II.  THIN  FILM  OXIDE  HETEROSTRUCTURES 


FIG.  1.  (Color  online.)  (a)  We  consider  oxide  het¬ 
erostructures  grown  along  the  (111)  direction  of  the  form 
AB’Oa/ABOa/AB’Oa  for  the  perovskite  structure  and 
AaB’aOy/AaBaOr/AaBbOT  for  the  pyrochlore  structure. 
Both  AB’Os  and  AaB’aOy  are  assumed  to  be  non-magnetic, 
wide  band  gap  (band)  insulators.  The  shaded  region  con¬ 
sists  of  a  (111)  bilayer  of  the  metallic  ABO3  perovskite,  such 
as  LaNiOs,  or  in  the  case  of  pyrochlore  structures,  a  bilayer 
or  a  trilayer  of  a  material  such  as  YalraOr.  (b)  Shown  are 
the  locations  of  the  transition-metal  ions  (B)  for  the  ABO3 
structure,  (c)  The  ABO3  bilayer  system  forms  a  “buckled 
honeycomb”  lattice,  which  appears  as  a  honeycomb  lattice 
when  projected  to  the  plane  perpendicular  to  (111).  We  as¬ 
sume  that  the  relevant  orbital  degrees  of  freedom  are  the  Cg 
orbitals  of  the  transition-metal  ions  for  ABO3.  For  A2B2O7 
with  heavy  transition  metal  ions,  the  relevant  orbitals  come 
from  the  spin-orbit  split  t2g  manifold  in  the  class  of  materials 
we  study  in  this  paper. 

We  are  interested  in  a  heterostructure  similar  to  that 
shown  in  Fig.l,  where  the  film  is  grown  along  the  (111) 
direction.  Most  experimental  studies  of  thin  films  in  the 
perovskite  systems  ABO3  are  grown  along  the  (001)  (or 
symmetry  related)  direction  because  aligning  the  growth 
direction  along  a  cubic  axis  generally  favors  high  quality 
films.  On  the  other  hand,  growing  along  (111)  effectively 
changes  the  crystal  structure  in  a  single  layer  thin  film 
from  a  square  lattice  of  transition  metal  ions  for  (001)  to 
a  triangular  lattice  of  transition  metal  ions  for  (111),  as 
shown  in  Fig.l.  For  a  bilayer  grown  along  (001),  one  will 
have  two  square  lattices  of  transition  metal  ions  stacked 
directly  on  top  of  each  other,  but  for  growth  along  (111) 
a  second  “shifted”  triangular  lattice  will  sit  on  top  of 
the  original  one.  The  combination  of  the  two  shifted 
triangular  lattices  is  a  “buckled”  honeycomb  lattice.  In 
this  example,  growing  along  the  (111)  direction  allows 
one  to  effectively  “engineer”  the  lattice  of  the  thin  film, 
with  important  implications  for  the  band  structure  in  the 
weak  coupling  limit  and  the  magnetic  order  in  the  strong 


coupling  limit. 

Likewise,  the  (111)  growth  direction  for  the  A2B2O7 
pyrochlores  leads  to  alternating  planes  of  triangular  and 
kagome  lattices  for  the  transition  metal  ions.  [See 
Fig. 5(a).]  To  the  best  of  our  knowledge,  there  have  been 
no  published  experimental  results  on  bilayer  or  trilayer 
films  of  A2B2O7  grown  along  (111),  though  there  are  a 
number  of  systems  where  (111)  growth  of  ABO3  films  has 
been  demonstrated, and  also  of  AB2O4  spinels. 
Various  theoretical  proposals  now  exist  for  topological 
phases  in  (111)  grown  transition  metal  oxide  films, 
and  the  list  of  candidate  materials  continues  to  grow.  We 
believe  it  is  likely  that  experiment  will  indeed  find  evi¬ 
dence  of  topological  phases  in  this  class  of  systems.  Once 
a  single  example  is  found,  experiments  can  be  done  to  op¬ 
timize  material  choices  by  “perturbing”  around  this  ma¬ 
terial  with  different  isovalent  elements,  substrates,  etc, 
in  order  to  achieve  the  maximum  bulk  gap.  Our  cal¬ 
culations  suggest  that  the  Chern  insulator  phase  stands 
out  as  the  mostly  likely  topological  candidate  for  realistic 
conditions  in  thin  film  oxide  heterostructures. 

To  be  concrete,  we  will  focus  on  the  LaNi03 
perovskite^^^^"*  and  the  ¥2^207  pyrchlore  iridate.^^’^^ 
(We  note  that  an  interesting  theoretical  study  of 
Co-doped  LaNi03  (111)  grown  bilayers  suggest  that 
correlation-driven  odd  parity  superconductivity  may  ap¬ 
pear  in  this  system.^®)  First  principles  calculations 
show  that  the  bands  around  the  Fermi  energy  are  pre¬ 
dominately  of  d-orbital  character  in  both  LaNi03  and 
Y2lr207.^^’^"*^’‘*®  As  a  result,  the  spatial  shape  of  the  or¬ 
bitals  are  highly  asymmetric  and  can  lead  to  interesting 
band  features  on  the  triangular,  kagome,  and  honeycomb 
lattices  that  appear  in  the  thin  film  structures  of  inter¬ 
est  to  us.  In  particular,  a  simple  tight-binding  model  for 
the  Ni  d-orbitals  leads  to  flat  bands  that  touch  dispersing 
bands  quadratically^^  (the  flat  touching  feature  persists 
with  more  realistic  first  principles  band  calculations^^’^^), 
and  also  bands  that  cross  in  a  Dirac  point  (see  Fig. 2). 
The  stability  of  such  band  touching  and  crossing  points 
with  respect  to  interactions  has  been  discussed  recently 
in  the  literature  in  the  context  of  interaction-generated 
topological  phases. 

The  central  idea  is  that  certain  types  of  interactions, 
such  as  an  on-site  multi-orbitaP^’^^  interaction  or  differ¬ 
ent  site  density-density  interaction  in  a  single  band^^^^'^ 
model,  can  generate  an  effective  spin-orbit  term  (at  the 
Hartree-Fock  level)  that  favors  topological  phases.  This 
is  one  of  the  ideas  we  will  explore  in  the  remainder  of  this 
article  in  the  context  of  real  materials,  though  we  will 
find  that  for  heavier  transition  elements,  such  as  irid¬ 
ium,  “interaction  generated  spin-orbit  coupling  terms” 
are  not  needed  to  access  topological  phases.  Thus,  any 
doubts  about  the  reliability  of  the  Hatree-Fock  predic¬ 
tions  for  enhanced  spin-orbit  effects  can  be  circumvented 
by  focusing  on  classes  of  materials  with  intrinsically  large 
spin-orbit  coupling. 
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A.  LaNiOs  bilayers 


FIG.  2.  First  principles  band  structure  of  the  fully  lattice  re¬ 
laxed  (LaA103)io/(LaNi03)2/(LaA103)io  system  within  (a) 
LDA  and  (b)  GGA,  which  are  essentially  indistinguishable. 
Note  the  quadratic  band  touching  at  the  F  point  as  well  as 
the  linear  crossings  at  K  and  For  a  spin  unpolarized  sys¬ 
tem  the  Fermi  energy  lies  right  at  the  quadratic  band  touching 
point,  while  for  a  fully  polarized  (ferromagnetic)  system  the 
Fermi  energy  lies  right  at  the  Dirac  point.  Our  Hartree-Fock 
calculations  suggest  that  for  realistic  interaction  values,  the 
system  is  very  close  to  a  fully  spin-polarized  state  with  a  quan¬ 
tized  Ghern  number-a  quantum  anomalous  Hall  state. 

See  Fig.4(b). 


The  band  structure  obtained  with  the  local  density 
approximation  (LDA)  and  the  generalized  gradient  ap¬ 
proximation  (GGA)^^  for  the  fully  relaxed  LaNiOs  (111) 
bilayer  is  shown  in  Fig.  2.  The  two  are  nearly  indis¬ 
tinguishable.  Rotations  of  the  octahedral  oxygen  cages 
are  known  to  be  important  for  large  classes  of  transi¬ 
tion  metal  oxides,®^  including  LaNiOs  for  which  the  adja¬ 
cent  oxygen  cages  counter-rotate  about  the  (111)  axis,  as 
shown  in  Fig. 3(a).  To  perform  Hartree-Fock  calculations 
with  this  band  structure,  we  consider  a  tight-binding 
model  based  only  on  the  nickel  Cg  orbitals  that  includes 
nearest-neighbor  hopping  via  the  oxygen  p-orbitals  and 
second-neighbor  hopping  via  the  oxygen  p-orbitals.  We 
find  a  a  good  fit  by  including  the  small  differences  in 
the  hopping  to  “outer”  versus  “inner”  oxygen  atoms. 
Assuming  trigonal  symmetry  is  preserved  (a  result  con¬ 
sistent  with  our  fully  relaxed  DFT  results),  we  take  the 
nearest-neighbor  Slater-Koster  parameters  for  hopping 
along  the  z-direction  to  be  described  by  the  matrix 


FIG.  3.  (Golor  online.)  (a)  The  pattern  of  the  octahe¬ 
dral  tilts/counter  rotations  present  in  the  fully  relaxed  struc¬ 
ture.  (b)  Layer  resolved  octahedral  rotation  angles  for  the 
(LaNi03)2/(LaA103)io  supercell  obtained  within  the  LDA 
approximation  to  DFT.  As  Fig.  2  shows,  these  rotations  do 
not  lift  the  quadratic  band  touching  at  the  F  point  or  the 
Dirac  points  at  K,  K’  in  the  Brillouin  zone  because  they  pre¬ 
serve  the  trigonal  point  group  symmetry.^"*  This,  in  turn,  im¬ 
plies  the  predictions  for  interaction-driven  topological  phases 
in  the  (LaNi03)2/(LaA103)jv  system  remain  qualitatively  un¬ 
changed  compared  to  the  “ideal”  lattice  structure.^® 


by  27r/3  is  .R  =  ^  ^72)  ’  ^  result,  we  find 

G  =  R^tzR,  ty  =  IRtxR.  The  Slater-Koster  parame¬ 
ters  for  second-neighbor  hopping  define  the  matrix^^ 


/  t'l2  V3A/2\ 
-‘it!  12  )  ■ 


(2) 


The  off-diagonal  entries  proportional  to  A  are  allowed  in 
the  bilayer  system  discussed  here  (as  opposed  to  a  perfect 
cubic  system)  because  the  two  possible  paths  connect¬ 
ing  second-neighbor  transition-metal  ions  are  not  equiv¬ 
alent:  they  either  involve  “inner”  or  “outer”  oxygens. 
Note  that  t^y  is  not  symmetric  if  A  yf  0  which  means 
that  there  is  an  associated  direction  for  the  hopping. 
We  use  the  convention  that  t^y  denotes  the  hopping 
of  an  electron  along  a  second  neighbor  bond  which  is 
reached  by  first  following  the  y-axis  and  then  the  cc-axis 
of  the  cube.  By  rotating  the  orbitals,  we  also  obtain 
the  second-neighbor  hopping  along  the  other  directions: 
tyz  R  txyR^  tzx  R  tyzR- 

The  generalized  tight-binding  model  now  takes  the 
form 


in  the  basis  {dz^ ,  dx^-y^).  Here  t  includes  predomi¬ 
nantly  the  hopping  via  the  intermediate  oxygen  while 
ts  arises  from  the  direct  overlap  and  is  small.  We  set 
fa  =  0  in  the  following.  Assuming  that  the  nearest- 
neighbor  hopping  in  the  x  and  y  directions  are  equiv¬ 
alent  to  the  hopping  along  the  z  direction,  we  obtain  the 
corresponding  matrices  by  a  rotation  of  the  Cg-orbitals 
around  (111)  by  ±27r/3.  The  matrix  for  the  rotation 


Hq  —  ^  ^  ^  ^  ^  ^  ^(^s^r^udsjV+eu  “f  h.C.^ 
r^A  s  u—xyz 

+  Y.Y.  S  (4  T*^u,Li+l^s,r+e^j,  — e^+i  ll.C.  )(3) 

r^A  s  U—xyz 

^  ^  ^  ^  ^  ^  (^^s,r^u,u-\-lds,r—eu-\-eu+i  “1“  ^L.C.^  , 

rG  B  s  U—xyz 

where  dg  =  {dz^^g,  dx2-y2  g)’^  is  a  vector  in  orbital  space, 
s  =t,  4-  is  the  spin  and  the  notation  u  -|-  1  refers  to  y  if 
u  =  X  with  a  cyclic  extension  to  the  other  elements. 
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fit  1 

t[eV] 

F  [eV] 

A  [eV] 

Ef  [eV] 

unrelaxed  (LDA) 

0.598 

0.062 

-0.023 

-0.693 

fully  relaxed  (LDA) 

0.541 

0.045 

-0.017 

-0.641 

fully  relaxed  (GGA)  | 

0.508 

0.046 

-0.016 

-0.593 

TABLE  I.  Parameters  obtained  in  tight-binding  fits  to  the  Cg 
DFT  band  structure  of  the  unrelaxed  and  fully  relaxed  super¬ 
lattice.  There  are  very  little  changes  between  the  unrelaxed 
and  relaxed  parameters,  with  t' /t  nearly  invariant  implying 
negligible  change  in  the  phase  diagram  shown  in  Fig. 4. 


Using  the  tight-binding  model  Hq  with  parameters  t, 
t' ,  and  A  (with  ts  =  0),  we  fitted  both  the  LDA  and 
GGA  band  structures  of  the  fully  relaxed  system  near  the 
Fermi  level.  The  fitting  parameters  are  listed  in  Table  I 
and  Fig.  4(a)  shows  the  LDA  together  with  the  tight- 
binding  band  structure  for  the  best  fit.  To  study  the 
multi-orbital  interactions  within  Hartree-Fock  theory,  we 
use  an  on-site  interaction  of  the  form^°’^^ 


Hint  = 


^rcni'^raj'  '^)  E  ^ras^rjSs 

r  a  a>0,s 


a//3  a//3 

a//3 


(4) 


We  assume  U'  =  U  —  2J  and  I  =  J,  which  are  valid 
in  free  space  and  believed  to  be  approximately  true  in 
the  solid  state  environment.  The  total  multi-orbital 
Hubbard  Hamiltonian  for  the  eg  electrons  is  given  by 
H  =  Hq  +  Ffint,  where  Hq  is  given  in  Eq.  (3).  The 
results^^’^"*^  are  shown  in  Fig. 4(b).  Studying  a  model 
that  explicitly  includes  oxygen  p-orbitals  (where  charge- 
transfer  physics  may  appear)  leads  to  similar  results. 


(a)  LDA  Tight  Binding  Fit 
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(b)  Hartree-Fock  Phase 
Diagram 


FIG.  4.  (color  online)  (a)  Fully  relaxed  LDA  band  structure 
and  tight-binding  (TB)  fit.  (b)  Hartree-Fock  phase  diagram 
of  the  LaNiOs  bilayer.  We  estimate  the  experimental  system 
has  parameters  t' /t  ~  0.1  and  J/U  ~  0.1  — 0.2,  which  is  rather 
close  to  the  quantum  anomalous  Hall  (QAH)  phase  with  ferro¬ 
magnetic  order.  FO=ferro-orbital,  AFM=antiferromagnetic, 
AFO=antiferro-orbital,  FM=ferromagnetic,  CDW=charge 
density  wave. 


B.  Y2lr207  bilayers  and  trilayers 


(a)  Pyrochlore  structure 


FIG.  5.  (Golor  online.)  (a)  Bulk  pyrochlore  lattice  structure 
showing  the  alternation  of  kagome  planes  (green  balls  on  lat¬ 
tice  sites)  and  triangular  lattice  planes  (grey  balls  on  lattice 
sites)  along  the  (111)  direction,  (b)  A  bilayer  film  viewed 
from  the  (111)  direction.  We  focus  on  the  bilayer  and  the 
triangular-kagome-triangular  structure,  which  show  the  most 
promise  for  topological  phases. 


In  order  to  study  the  (111)  films  of  the  pyrochlore  in 
a  set-up  similar  to  that  shown  in  Fig.  1(a),  we  consider 
a  tight-binding  model 

Hq  =  ^  ^  ^iaJpCic^Cjp  —  A  ^  (  1;  •  Sj,  (5) 

{id),a,P  * 

where  the  d-orbital  hopping  takes  the  form^®’®^ 

tiaJP  =  tlaJiS  +  (6) 

which  contains  both  an  indirect  and  a  direct  hopping 
term  between  the  d-orbitals. ^^43,55  A  >  0  is  the 

intrinsic  spin-orbit  coupling  in  the  system  which  acts 
within  the  t2g  manifold  so  |1|  =  1,  and  is  the  spin  of  the 
electron  in  a  t2g  d-orbital  on  site  jA^43.55 
ides,  the  strength  of  the  spin-orbit  coupling  is  estimated 
to  be  0.2-0. 7  eV  and  the  hopping  strength  is  on  the  or¬ 
der  of  0.4-0. 6  eV.^^40  q^j^g  hopping  amplitude  in  Eq.  (6) 
contains  a  direct  d-d  hopping,  in  addition  to  the 

indirect  hopping  via  the  oxygen  orbitals. xhe  direct 
hopping  is  parameterized  by  the  strength  of  the  cr-bonds, 
ts,  and  the  7r-bonds,  tp.  Following  Refs.  [19]  and  [54],  we 
consider  a  set  of  representative  ratios  to  explore  a  real¬ 
istic  parameter  space:  We  set  tp  =  —2ts/3  and  consider 
the  cases  of  =  —t  and  tg  =  t.  Our  preliminary  GGA 
calculations  for  the  thin-films  suggest  the  band  structure 
is  most  similar  to  that  for  tg  =  —t,  though  the  further 
neighbor  hopping  is  considerably  more  important  for  the 
5d  orbitals  in  ¥2^207  than  it  is  for  the  3d  orbitals  in 
LaNiOg. 

To  carry  out  Hartree-Fock  calculations  for  the  Y2lr207 
films,  we  use  the  tight-binding  model  in  Eq.(5)  with 
tg  =  —t  supplemented  by  an  on-site  Hubbard  term, 
Hjj  =  ^14*44,  and  further  restrict  ourselves  to  the 

j  =  1/2  manifold.^®  The  results  are  shown  in  Fig.  6.  Note 
the  triangular-kagome-triangular  (TKT)  film  supports  a 
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fairly  wide  region  of  a  quantum  anomalous  Hall  (Chern 
insulator=CI)  state.  We  find  that  if  one  includes  fluctu¬ 
ations  beyond  the  Hartree-Fock  approximation,  this  re¬ 
gion  moves  to  larger  U  values,  close  to  what  is  reason¬ 
able  for  Y2lr207.  Therefore,  we  conclude  that  both  the 
(111)  grown  LaNiOa  bilayer  and  the  (111)  grown  Y2lr207 
TKT  trilayer  are  candidates  for  a  quantum  anomalous 
Hall  state.  Moreover,  we  And  that  small  changes  to  the 
kinetic  terms  in  Eq.(3)  can  lead  to  a  Z2  topological  insu¬ 
lator  in  the  (111)  grown  Y2lr207  bilayer  for  small  U 


U/t 

FIG.  6.  (Color  online.)  Hartree-Fock  phase  diagrams  for  the 
single  kagome  layer  (monolayer),  bilayer,  triangular-kagome- 
triangular  (TKT),  and  kagome-triangular-kagome  (KTK)  sys¬ 
tems.  We  have  used  =  —t  and  restricted  ourselves 
to  the  j  =  1/2  manifold.  M=metallic,  I=trivial  insu¬ 
lator,  TI=topological  insulator,  MC=magnetic  conductor, 
MI=magnetic  insulator,  CI=Chern  insulator  (same  as  QAH). 


III.  CONCLUSIONS 

In  conclusion,  we  have  shown  both  the  (111)  grown 
LaNiOa  bilayer  and  the  (111)  grown  Y2lr207  TKT  tri¬ 
layer  are  candidates  for  a  quantum  anomalous  Hall  state, 
which  will  show  a  quantize  Hall  conductance.  We  also 
found  that  realistic  changes  to  the  kinetic  terms  of 
Y2lr207  can  also  lead  to  a  Z2  topological  insulator  in 
the  (111)  grown  Y2lr207  bilayer.  The  most  natural  way 
to  detect  these  states  is  through  transport  measurements, 
as  has  been  done  for  the  few  known  Z2  and  QAH 

systems.®®  From  the  point-of-view  of  technology  applica¬ 
tions,  there  is  clearly  a  need  in  the  held  for  more  systems 
known  to  possess  these  phases.  The  huge  variety  of  tran¬ 
sition  metal  oxides,  combined  with  existing  theoretical 
results,  would  seem  to  suggest  it  may  only  be  a  matter  of 
time  before  a  topological  phase  is  discovered  in  a  system 
similar  to  those  we  considered  here. 
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Topological  insulators  have  emerged  as  a  major  topic  of  condensed  matter  physics  research  with 
several  novel  applications  proposed.  Although  there  are  now  a  number  of  established  experimental 
examples  of  materials  in  this  class,  all  of  them  can  be  described  by  theories  based  on  electronic 
band  structure,  which  implies  that  they  do  not  possess  electronic  correlations  strong  enough  to 
fundamentally  change  this  theoretical  description.  Here,  we  review  recent  theoretical  progress  in 
the  description  of  a  class  of  strongly  correlated  topological  insulators  —  fractionalized  topological 
insulators  —  where  band  theory  fails  dramatically  due  to  the  fractionalization  of  the  electron  into 
other  degrees  of  freedom. 


The  field  of  topological  insulators  (TI)  is  now  approx¬ 
imately  a  decade  old,  and  there  are  a  number  of  excel¬ 
lent  reviews  available,  from  both  an  experimental  and 
theoretical  perspective^^^.  Our  goal  in  this  article  is  to 
provide  a  progress  report  on  a  specihc  subdirection  that 
has  developed  rapidly  in  the  last  few  years — the  theo¬ 
retical  study  of  phases  of  matter  that,  while  conceptu¬ 
ally  related  to  topological  insulators,  cannot  be  described 
by  an  electronic  band  structure  because  of  the  fraction¬ 
alization  of  the  electron  into  other  degrees  of  freedom, 
driven  by  strong  electronic  correlations.  This  latter  sce¬ 
nario  has  earlier  precedents  in  the  Mott  insulator  (  which 
exhibits  the  separation  of  spin  and  charge)  and  the  frac¬ 
tional  quantum  Hall  effect  (which  has  excitations  with 
only  a  fraction  of  the  charge  of  the  electron),  but  ac¬ 
quires  a  distinct  flavor  when  applied  to  topological  in¬ 
sulators.  By  taking  fractionalized  topological  insulators 
as  an  example,  we  hope  to  convey  the  broader  message 
that  the  class  of  materials  commonly  known  as  topologi¬ 
cal  insulators  barely  scratches  the  surface  of  the  possible 
types  of  topological  materials,  once  electron  correlations 
are  taken  into  account. 


SPLITTING  THE  ELECTRON  APART 

We  restrict  our  discussion  to  two  and  three  spatial  di¬ 
mensions,  as  fractionalized  excitations  in  gapped  one¬ 
dimensional  systems  are  require  additional  considera¬ 
tions.  Since  we  focus  on  systems  not  described  by  band 
theory,  we  merely  note  in  passing  that  interactions  can 
drive  a  topologically  trivial  phase  into  a  topological  band 
insulator  through  the  spontaneous  generation  of  spin- 
orbit  coupling®’®  or  magnetic  order^,  or  through  the 
renormalization  of  electronic  bands®.  By  contrast,  in 
the  phases  we  highlight  below  the  quantum  numbers  of 
the  electron  break  apart,  implying  that  the  latter  is  not 
a  sharply  defined  quasiparticle  and  the  system  cannot 
be  adequately  described  by  an  electronic  band  structure. 
The  chief  theoretical  challenges  for  the  study  of  fraction¬ 
alized  topological  insulators  are:  (i)  to  describe  them  in 
a  conceptually  convenient  way,  (ii)  to  predict  the  condi¬ 


tions  under  which  they  are  expected  (i.e.,  find  the  set  of 
Hamiltonians  that  support  them) ,  and  if  possible  (iii)  to 
provide  a  list  of  specific  candidate  materials  that  can  be 
investigated  experimentally.  The  field  is  presently  at  the 
stage  of  tackling  challenges  (i)  and  (ii),  with  relatively 
little  work  on  (iii). 

A  convenient  way  to  describe  the  basic  properties  of 
fractionalized  topological  insulators  is  the  slave-particle 
approach.  In  the  present  context,  it  answers  the  follow¬ 
ing  question:  What  new  state  of  matter  does  one  get  if 
electrons  inside  a  topological  insulator  fractionalize  into 
other  degrees  of  freedom?  The  answer  depends  on  the 
type  of  topological  insulator  with  which  one  starts,  as 
well  as  the  assumed  pattern  of  electron  fractionalization. 
From  these  basic  inputs,  the  slave-particle  approach  pre¬ 
dicts  the  universal  properties  of  the  resulting  fraction¬ 
alized  topological  insulator,  such  as  the  quantum  num¬ 
bers  and  quantum  statistics  of  low-energy  excitations, 
the  degeneracy  of  the  ground  state  on  surfaces  of  vari¬ 
ous  topologies,  and  quantized  electromagnetic,  spin,  and 
thermal  response  properties  (Table  I) .  The  slave-particle 
approach  does  not,  however,  solve  the  much  more  diffi¬ 
cult  problem  of  the  classification  of  all  possible  correlated 
topological  phases — which  goes  beyond  the  scope  of  our 
focused  review — because  a  generic  correlated  topological 
phase  may  not  always  be  understood  as  originating  via 
fractionalization  from  a  noninteracting  topological  insu¬ 
lator.  More  sophisticated  methods  such  as  the  group 
cohomology®’^®  and  cobordism^^  approaches  have  been 
recently  developed  to  address  this  considerable  challenge. 

Besides  fractionalization,  to  obtain  an  exotic  topologi¬ 
cal  phase  one  typically  assumes  that  some  of  the  fraction¬ 
alized  degrees  of  freedom  form  a  topological  band  insu¬ 
lator.  The  latter  can  thus  be  seen  as  a  nonfractionalized 
parent  state  from  which  fractionalized  daughter  states 
can  be  constructed  (Table  I,  first  two  columns).  The 
nontrivial  response  properties  of  the  parent  state  combine 
with  the  fractional  quantum  numbers  of  the  fractional¬ 
ized  degrees  of  freedom  to  yield  even  more  exotic  response 
properties  for  the  daughter  state.  This  scenario  occurs 
in  slave-particle  mean-field  studies  of  models  of  topologi¬ 
cal  band  insulators  subjected  to  strong  electron-electron 
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parent  state 

fractionalized  state 

d 

symmetry 

EM  response 

spin  response 

thermal  response 

Cl 

CSL 

2 

^xy  ~  0 

=  0 

KxylT  —  2 

QSH 

fractionalized  QSH 

2 

TRS 

^xy  —  0 

G%y  =  0 

Hxy/T  =  0 

TI/WTI 

TMI/WTMI 

3 

TRS 

(9  =  0 

Sgrav  =  TT  (TMI) 

TCI 

TCMI 

3 

mirror 

Cl 

Cl* 

2 

^xy  —  2 

al,y  =  0 

KxyjT  —  2 

QSH 

QSH* 

2 

TRS 

^ xy  - —  0 

S  _  C\ 

^xy  ^ 

HxyjT  =  0 

TI 

TI* 

3 

TRS 

6  =  TV 

^grav  —  rr 

Cl 

FCI 

2 

^ xy  ~  1/3 

G%y  =  0 

K-xyjT  =  1 

QSH 

FQSH 

2 

TRS 

^  xy  —  0 

o%y  =  2/3 

Hxy/T  =  0 

TI 

FTI 

3 

TRS 

Q  =  7r/3 

^grav  —  TT 

TABLE  I.  Summary  of  the  d-dimensional  fractionalized  topological  insulators  discussed  in  this  article  and  their  universal 
response  properties  (acronyms  are  defined  in  the  main  text),  a^y'-  Hall  conductance  in  units  of  ^ lh\  9:  coefficient  of  the  E  ■  B 
term  in  units  of  e^/27rh;  a%y-.  spin  Hall  conductance  in  units  of  e/47r  for  a  system  with  conserved  2;  component  of  spin;  Kxy/T-. 
thermal  Hall  conductance  divided  by  the  temperature,  in  units  of  dgrav^  coefficient  of  the  gravitational  equivalent  of 

the  E  ■  B  term,  responsible  for  a  surface  quantized  thermal  Hall  effect. 


interactions.  Besides  interacting  with  external  fields  via 
their  fractional  quantum  numbers  such  as  spin  or  charge, 
the  fractionalized  degrees  of  freedom  also  interact  with 
emergent  gauge  fields.  As  in  quantum  chromodynamics 
(QCD),  these  gauge  fields  are  confining  in  a  conventional 
phase  where  fractionalized  degrees  of  freedom — “quarks” 
in  the  QCD  analogy — do  not  exist  at  low  energies  as 


FRACTIONALIZATION  AND  NEW  PHASES  OF 
MATTER 

The  hrst  pattern  of  electron  fractionalization  to  be  con¬ 
sidered  in  the  context  of  strongly  correlated  topological 
insulators  was  spin-charge  separation,  originated  in  the 
study  of  the  Luttinger  liquid  in  ID  and  the  Mott  tran¬ 
sition  in  2D^^d3  electron  is  assumed  to  fractional¬ 

ize  into  a  charged,  spinless,  bosonic  rotor  and  an  elec¬ 
trically  neutral,  spinful,  fermionic  spinon,  coupled  by  a 
U{1)  gauge  field  (Fig.  1).  The  electronic  quasiparticle 
weight  is  proportional  to  the  expectation  value  of  the 
rotor  field.  The  resulting  theory  typically  has  two  dis¬ 
tinct  stable  phases:  a  weakly  interacting  phase  where 
the  rotors  undergo  Bose  condensation  and  the  rotor  field 
acquires  a  nonzero  expectation  value,  and  a  strongly  in¬ 
teracting  phase  where  the  rotors  are  uncondensed.  In 
the  weakly  interacting  phase,  the  electronic  quasiparti¬ 
cle  weight  is  hnite  and  the  electronic  band  structure  is 
well  defined,  while  in  the  strongly  interacting  phase  the 
quasiparticle  weight  vanishes  due  to  strong  rotor  fluc¬ 
tuations,  corresponding  to  a  Mott  phase.  In  the  latter 
case,  the  spinons  acquire  a  band  structure  of  their  own 
which,  if  gapped,  can  be  classified  by  the  usual  topologi¬ 
cal  invariants'’^.  In  2D,  Chern-insulator  (Cl)  or  Z2  topo¬ 
logical  insulator  band  structures  for  the  spinons  lead  to 
a  chiral  spin  liquid^^  (CSL)  or  a  fractionalized  quantum 


long-lived,  propagating  excitations.  Inside  a  fractional¬ 
ized  topological  insulator,  the  gauge  fields  are  deconfined 
and  allow  the  fractionalized  degrees  of  freedom  to  lead  an 
existence  of  their  own.  The  gauge  fields  themselves  ad¬ 
mit  propagating  excitations — the  gauge  bosons — which 
can  be  gapless  or  gapped,  depending  on  the  type  of  gauge 
field  and  the  nature  of  the  parent  topological  phase. 


spin  Hall  (QSH)  insulator^®’^®,  respectively.  The  frac¬ 
tionalized  QSH  insulator  is  only  stable  in  the  presence 
of  additional  layers  containing  enough  gapless  spinons  to 
screen  the  U{1)  gauge  field  and  avoid  the  instanton  effect 
that  causes  confinement^®.  In  3D  time-reversal  invari¬ 
ant  systems,  the  Mott  phases  are  called  strong  topolog¬ 
ical  Mott  insulators^^  (TMI)  or  weak  topological  Mott 
insulators^®  (WTMI)  depending  on  whether  the  spinons 
form  a  strong  or  weak  Z2  topological  insulator,  respec¬ 
tively.  The  same  analysis  applied  to  systems  with  mirror 
crystalline  symmetries  leads  to  the  concept  of  a  topo¬ 
logical  crystalline  Mott  insulator^®  (TCMI).  All  these 
3D  fractionalized  topological  phases  have  gapless  sur¬ 
face  modes  of  spinons  with  the  properties  that  electrons 
would  have  in  the  corresponding  noninteracting  topolog¬ 
ical  phases,  but  in  addition  have  a  bulk  gapless  “photon” 
mode  from  the  emergent  U{1)  gauge  field^^’^°.  How¬ 
ever,  there  is  no  quantized  electromagnetic  response  of 
the  E  ■  B  type^“®  because  the  spinons  are  electrically 
neutral.  To  identify  this  class  of  spin-charge  separated 
states  experimentally,  methods  similar  to  those  proposed 
or  used  to  identify  quantum  spin  liquids^^  should  be  ap¬ 
plied.  One  should  observe  a  bulk  contribution  to  the  low- 
temperature  specific  heat  arising  from  the  gapless  photon 
mode,  although  it  could  be  difficult  to  distinguish  from 
the  phonon  contribution.  In  addition,  one  could  imagine 
detecting  the  linearly  dispersing  surface  spinon  modes  via 
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FIG.  1.  Fractionalized  topological  insulators  accessed  by 
the  U{1)  slave-rotor  description.  Outside  the  material,  the 
charged  bosonic  rotor  (blue  sphere)  and  the  neutral  fermionic 
spinon  (red  sphere)  are  confined  together  to  form  an  elec¬ 
tron.  Inside  the  material,  they  become  deconfined  particles 
and  interact  via  an  emergent  U(l)  gauge  field  (green  wiggly 
line).  In  the  2D  chiral  spin  liquid  (a),  these  particles  acquire 
semionic  (i.e.,  half-fermionic)  statistics  and  there  are  chiral 
edge  states  of  spinons.  In  the  2D  fractionalized  quantum  spin 
Hall  insulator  (b),  the  spinon  edge  states  are  helical.  In  the 
3D  topological  Mott  insulator  states  (c),  the  U{1)  gauge  field 
is  gapless  in  the  bulk,  and  there  are  gapless  Dirac-like  spinon 
surface  states. 

the  observation  of  2D  spin  or  heat  transport  localized  at 
the  surface  of  the  material. 

A  different  pattern  of  electron  fractionaliza- 
tion  originating  from  theories  of  high-temperature 
superconductivity^^  assumes  that  the  electron  fraction¬ 
alizes  into  a  neutral,  spinless,  bosonic  slave-Ising  particle 
and  a  charged,  spinful  slave-fermion,  coupled  by  a  Z2 
gauge  field^^  (Fig-  2).  This  type  of  fractionalization  may 
appear  trivial,  as  the  slave-Ising  particle  carries  neither 
the  charge  nor  the  spin  of  the  electron.  The  resulting 
fractionalized  phases,  however,  are  far  from  trivial. 
As  in  the  C/(l)  theory,  the  Z2  theory  typically  has  a 
non-fractionalized  weakly  interacting  phase  where  the 
slave-Ising  particles  Bose  condense,  and  a  fractionalized 
strongly  interacting  phase  where  they  are  uncondensed. 
By  contrast  with  the  U (1)  theory  however,  in  the  Z2  case 
the  slave-fermion  carries  the  charge  of  the  electron,  thus 
the  fractionalized  phases  inherit  the  electromagnetic 
response  properties  of  the  parent  topological  band 
insulator.  In  2D,  Chern-insulator  and  Z2  topological 
insulator  band  structures  for  the  slave  fermions  lead  to  a 
correlated  Chern  insulator^'^  (Cl*)  with  quantized  Hall 
response  and  a  correlated  quantum  spin  Hall  insulator^® 
(QSH*),  respectively.  These  are  fully  gapped  phases 
with  intrinsic  topological  order  of  the  double-semion 
and  toric-code  type,  respectively.  In  3D  time-reversal 
invariant  systems,  a  Z2  topological  insulator  band 
structure  leads  to  a  strongly  correlated  topological 
insulator^®  (TI*)  that,  unlike  the  TMI  described  above, 
exhibits  a  quantized  E  ■  B  electromagnetic  response 
and  has  no  gapless  bulk  excitations.  The  TI*  is  a  true 


FIG.  2.  Fractionalized  topological  insulators  accessed  by  the 
Z2  slave-Ising  description.  Outside  the  material,  the  neutral 
bosonic  slave-Ising  particle  (orange  sphere)  and  the  charged 
slave-fermion  (purple  sphere)  are  confined  together  to  form 
an  electron.  Inside  the  material,  they  become  deconfined  par¬ 
ticles  and  interact  via  an  emergent  Z2  gauge  field  (orange 
wiggly  line)  that  supports  yr-flux  vortex  excitations  (orange 
disk).  In  the  2D  correlated  Ghern  insulator  (a),  vortices  ac¬ 
quire  semionic  statistics  and  there  are  chiral  edge  states  of 
slave-fermions.  In  the  2D  correlated  quantum  spin  Hall  insu¬ 
lator  (b),  vortices  have  semionic  mutual  statistics  with  slave- 
Ising  particles  and  slave- fermions,  and  the  slave- fermion  edge 
states  are  helical.  In  the  3D  correlated  topological  insulator 
(c),  the  TT-flux  vortex  becomes  a  bulk  vortex  loop,  and  there 
are  gapless  Dirac-like  slave-fermion  surface  states. 

3D  topological  phase  whose  bosonic  sector  is  described 
by  a  topological  BF  theory^^’^®  that  imparts  a  relative 
statistical  angle  of  tt  between  bulk  particle  and  loop 
excitations  and  is  responsible  for  a  nontrivial  topological 
ground-state  degeneracy. 

Finally,  there  are  even  more  exotic  fractionalized 
phases  in  which  the  electron  with  charge  e  splits  apart 
into  three  fermionic  partons  of  charge  e/3  (Fig.  3).  In  2D, 
when  these  partons  form  noninteracting  integer  quantum 
Hall  states  one  obtains  the  familiar  fractional  quantum 
Hall  states.  In  the  conceptually  similar  case  of  a  Chern 
insulating  band  structure,  one  obtains  fractional  Chern 
insulators  (FCI) — fractional  quantum  Hall  states  that  ex¬ 
ist  in  the  absence  of  an  external  magnetic  field.  There 
is  already  a  sizeable  literature  on  FCI,  and  we  direct  the 
reader  to  existing  review  articles  on  the  subject^®’®®  for 
more  information  and  original  references.  If  the  frac¬ 
tionally  charged  fermions  form  a  time-reversal  invari¬ 
ant  Z2  topological  insulator,  one  obtains  a  fractional 
quantum  spin  Hall  insulator  or  2D  fractional  topologi¬ 
cal  insulator®!-®®  (FTI),  and  in  3D  a  3D  FTI®®-®®.  The 
relevant  gauge  theories  in  the  example  above  are  of  the 
U{1)  X  [/(I)  or  Z3  variety,  leading  to  gapless  bulk  pho¬ 
tons  or  no  bulk  gapless  excitations,  respectively  (a  non- 
Abelian  SU{3)  gauge  theory  is  formally  possible  but  is 
typically  plagued  by  the  problem  of  confinement).  In 
both  cases  the  3D  FTI  exhibits  a  quantized  E  ■  B  electro¬ 
magnetic  response,  but  with  a  coefficient  one-third  that 
of  a  topological  band  insulator.  From  an  experimental 
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FIG.  3.  Fractional  topological  insulators.  Outside  the  ma¬ 
terial,  the  three  fermionic  partons  of  charge  e/3  are  confined 
together  to  form  an  electron.  Inside  the  material  and  in  the 
17(1)  X  17(1)  description,  they  become  deconfined  particles 
and  interact  via  two  17(1)  gauge  fields  (orange  and  magenta 
wiggly  lines).  In  the  2D  phases  (a)  and  (b),  the  partons  ac¬ 
quire  fractional  statistics.  In  the  fractional  Chern  insulator 
(a),  there  are  chiral  edge  states  of  partons,  while  in  the  frac¬ 
tional  quantum  spin  Hall  Insulator  (b)  the  parton  edge  states 
are  helical.  In  the  3D  fractional  topological  insulator  (c),  the 
1/(1)  gauge  fields  are  gapless  in  the  bulk,  and  there  are  gap¬ 
less  Dirac-like  parton  surface  states.  In  the  Z3  description, 
the  gauge  held  would  be  gapped  in  the  bulk  with  vortex  loop 
excitations  similar  to  those  of  Fig.  2(c). 

standpoint,  a  3D  FTI  should  be  easier  to  detect  than  a 
3D  TMI  because  its  gapless  surface  states  are  electrically 
charged,  leading  to  nontrivial  signatures  in  photoemis¬ 
sion,  tunneling,  and  quantum  oscillations^'^. 

Can  the  fractionalized  topological  insulators  whose 
phenomenology  we  have  described  so  far  be  reliably  ob¬ 
tained  as  ground  states  of  interacting  electronic  Hamil¬ 


tonians?  The  answer  is  a  partial  yes.  In  2D,  FCI  ground 
states  have  been  reliably  obtained  as  ground  states  of 
models  of  spinless  electrons  with  nearest-neighbor  repul¬ 
sive  interactions  via  numerical  studies^®’^°  and  strong¬ 
coupling  expansions'^^.  2D  time-reversal  invariant  FTI 
ground  states  have  likewise  been  found  in  numerical 
studies'^^’^'^  and  exactly  soluble  models^^’^®  of  interacting 
fermions  and/or  bosons  (compressible,  bosonic  analogs  of 
2D  FTI  may  have  also  been  found  in  microscopic  models 
of  interacting  bosons"^®).  Progress  has  also  been  made  in 
systematically  designing  Hamiltonians  for  which  various 
spin-charge  separated  states  such  as  the  CSL  are  exact 
ground  states^^’''®,  but  these  are  models  of  interacting  lo¬ 
calized  spins  rather  than  models  of  interacting  electrons. 
The  challenge  remains  mostly  open  in  3D — where  numer¬ 
ical  studies  come  at  a  prohibitive  computational  cost — 
with  the  exception  of  exactly  soluble  models^"*. 

We  have  discussed  a  few  examples  of  fractionalized 
topological  insulators.  While  a  list  of  specific  mate¬ 
rials  that  would  harbor  such  phases  is  currently  lack¬ 
ing,  materials  that  combine  the  fundamental  ingre¬ 
dients  of  strong  spin-orbit  coupling  and  strong  elec¬ 
tron  correlations — such  as  the  pyrochlore  iridates^®  and 
samarium  hexaboride®*^ — do  exist  and  are  currently  un¬ 
der  intense  experimental  study.  Such  materials,  as  well  as 
systems  of  ultracold  atoms  subjected  to  synthetic  gauge 
fields,  may  unveil  entirely  novel  quantum  phases  of  mat¬ 
ter,  topological  or  otherwise. 
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In  this  paper  we  review  some  connections  recently  discovered  between  topological  insulators  and 
certain  classes  of  quantum  spin  liquids,  focusing  on  two  and  three  spatial  dimensions.  In  two 
dimensions  we  show  the  integer  quantum  Hall  effect  plays  a  key  role  in  relating  topological  insulators 
and  chiral  spin  liquids  described  by  fermionic  excitations,  and  we  describe  a  procedure  for  “generating" 
a  certain  class  of  topological  states.  In  three  dimensions  we  discuss  interesting  relationships  between 
certain  quantum  spin  liquids  and  interacting  “exotic"  variants  of  topological  insulators.  We  focus 
attention  on  better  understanding  interactions  in  topological  insulators,  and  the  phases  nearby  in 
parameter  space  that  might  result  from  moderate  to  strong  interactions  in  the  presence  of  strong  spin- 
orbit  coupling.  We  stress  that  oxides  with  heavy  transition  metal  ions,  which  often  host  a  competition 
between  electron  interactions  and  spin-orbit  coupling,  are  an  excellent  place  to  search  for  unusual 
topological  phenomena  and  other  unconventional  phases, 
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1.  Introduction 

At  first  glance  topological  Insulators  and  quantum  spin  liquids 
have  nothing  to  do  with  each  other.  A  topological  insulator  (Tl)  has 
properties  that  can  be  understood  in  the  absence  of  any  electron 
Interactions  [1-3],  while  quantum  spin  liquids  (QSL)  are  largely 
understood  theoretically  In  terms  of  lattice  models  of  Interacting 
spins  [4],  Moreover,  topological  Insulators  explicitly  possess  charge 
and  spin  degrees  of  freedom,  while  spin  liquids  are  understood  only 
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in  terms  of  spin  degrees  of  freedom.  Quantum  spin  liquids  come  in 
both  gapped  and  gapless  varieties,  but  even  for  the  gapped  variety 
with  topological  order  it  would  appear  there  is  little  to  say  by  way  of 
a  relationship  between  TIs  and  QSLs,  As  a  supporting  argument,  one 
might  point  out  that  the  “type"  of  topological  order  is  different: 
Topological  insulators  do  not  have  any  non-trivial  ground  state 
degeneracy,  while  some  of  the  best  understood  topological  spin 
liquids  possess  a  non-trivial  ground  state  degeneracy  [5], 

However,  to  dismiss  any  relationship  between  these  two 
phases  of  matter  is  to  ignore  some  beautiful,  and  we  believe 
insightful,  relationships.  To  help  motivate  why  it  is  useful  to  look 
for  connections  between  systems  that  apparently  have  no 
relationship  to  one  another  we  need  to  look  no  further  than 
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Laughlin’s  account  of  his  discovery  of  the  fractionai  quantum  Hail 
effect  [6].  According  to  Laughiin,  a  one-dimensional  polymer  known 
as  polyacetylene  (which  displays  charge  fractionalization,  at  least 
theoretically  [7])  was  an  important  motivation  for  his  inspired  guess 
for  a  fractionai  quantum  Hall  state  wavefunction  [6].  As  it  so  turns 
out,  there  are  also  deep  connections  between  polyacetylene  and 
topological  insulators  [8,9],  We  believe  analogous  relationships 
between  TIs  and  QSLs  can  also  provide  useful  insights,  and  could 
play  a  key  role  in  deepening  our  understanding  of  strongly  inter¬ 
acting  variants  of  topological  insulators  that  are  not  adiabatically 
connected  to  the  non-interacting  Tl  limit. 

In  this  topical  review  we  will  flush  out  some  of  the  connections 
between  TIs  and  QSLs  recently  established.  In  this  paper,  we  will 
focus  on  a  number  of  different  lattice  models  of  fermions.  Some  of 
the  models  are  non-interacting,  some  are  interacting,  some  have  no 
intrinsic  spin-orbit  coupling,  while  some  do  have  intrinsic  spin-orbit 
coupling.  We  study  a  variety  of  lattices  in  two  and  three  spatial 
dimensions  and  focus  on  the  interplay  of  lattice  geometry,  interac¬ 
tions,  and  spin-orbit  coupling.  Our  aim  is  to  enlarge  our  under¬ 
standing  of  the  phase  diagram  of  interacting  fermions  in  two  and 
three  dimensions  by  solving  different  model  problems  and  using 
them  as  points  of  generalization.  For  example,  which  features  appear 
often  and  therefore  are  rather  general,  and  which  features  appear  to 
be  model  specific?  This  approach  is  complementary  to  directly 
develop  an  effective  theory.  Effective  theories  are  extremely  power¬ 
ful  because  of  their  generality  and  efficient  exploitation  of  the 
various  symmetries  present  in  a  problem.  However,  it  is  not  always 
straightforward  to  “guess”  an  effective  theory  of  a  new  phase  of 
matter  so  special  cases  derived  from  more  microscopic  lattice 
models  can  serve  as  a  useful  guide,  particularly  when  strong 
interactions  and/or  correlations  are  present.  That  is  the  point-of- 
view  we  take  in  this  paper. 

A  summary  of  our  main  results  is  as  follows.  For  the  case  of  two- 
dimensional  lattices  we  show  that  under  rather  general  conditions 
there  is  a  connection  between  TIs  and  Kitaev-type  [10,11]  spin 
models.  An  earlier  work  by  Lee  et  al.  [12]  connecting  tbe  honeycomb 
Kitaev  model  and  the  integer  quantum  Hall  effect  was  an  important 
inspiration  for  our  studies.  A  relationship  between  the  integer 
quantum  Hall  effect,  TIs,  and  Kitaev  models  exists  because  the  latter 
can  be  exactly  solved  by  a  mapping  that  reduces  the  interacting  spin 
degrees  of  freedom  to  non-interacting  Majorana  degrees  of  freedom 
[11[.  These  Majorana  fermions  hop  on  the  lattice  and  have  a  band 
structure  similar  to  “normal”  fermions  (only  with  a  redundancy  in 
the  labeling  of  the  states).  Therefore,  one  may  study  the  topological 
properties  of  the  energy  bands,  such  as  whether  they  possess  a  finite 
Chern  number.  Energy  bands  of  finite  Chem  number  will  turn  out  to 
establish  a  connection  between  TIs  and  Kitaev-type  spin  liquids  (and 
of  course  tbe  integer  quantum  Hall  effect). 

Staying  with  two-dimensions,  we  will  discuss  recent  results  on 
obtaining  a  fractional  quantum  Hall  effect  in  a  fiat-band  lattice 
model  with  a  partially  filled  band  with  a  finite  Chern  number.  This 
lattice  model  may  prove  useful  in  the  studies  of  “fractional” 
topological  insulators  in  two-dimensions  with  time-reversal  sym¬ 
metry  [13-16].  We  also  describe  a  class  of  models  in  which  spin- 
orbit  coupling  can  be  spontaneously  generated  at  the  mean-field 
level  [17].  We  contrast  the  relative  instability  to  electron  interac¬ 
tions  of  Dirac-like  and  quadratic  band  touching  points  in  the  non¬ 
interacting  band  structure.  We  find  that  when  quadratic  band 
touching  points  are  present,  topological  phases  often  appear  as  the 
leading  instability  with  interactions.  Such  mean-field  studies  reduce 
a  many-particle  problem  to  a  single-particle  problem  and  serve  at 
least  two  key  roles  in  the  study  of  topological  insulators:  (1)  They 
elucidate  where  topological  phases  can  be  expected  to  appear  in  the 
presence  of  electron  interactions  and  (2)  they  help  guide  the  class  of 
non-interacting  lattice  models  that  contain  interesting  physics.  We 
finish  our  discussion  of  two-dimensional  systems  with  an  example 


of  theoretical  Tl  studies  guiding  the  discovery  of  a  QSL  with  novel 
properties  [18]  and  comment  on  the  role  that  disorder  can  play  in 
stabilizing  topological  quantum  spin  liquids  [19]  and  topological 
insulators  [20-24]. 

In  tbe  case  of  three-dimensions,  we  focus  on  systems  with 
intrinsic  spin-orbit  coupling  and  intermediate  strength  interactions. 
With  an  eye  towards  transition  metal  oxides  with  4d  and  5d 
elements,  we  use  the  slave-rotor  mean-field  theory  [25,26]  to 
investigate  the  competition  between  interactions,  spin-orbit  cou¬ 
pling  and  lattice  distortions  on  the  pyrochlore  lattice.  We  find  that  a 
new  phase — the  “weak  topological  Mott  insulator” — appears  [27] 
with  topologically  protected  gapless  modes  that  carry  heat  only 
along  a  certain  class  of  bulk  defects  ]28].  As  a  function  of  pressure,  a 
complex  set  of  transitions  between  conventional  and  exotic  phases 
is  possible  in  such  systems.  We  discuss  the  connection  of  the  slave- 
rotor  mean-field  states  to  quantum  spin  liquids  describe  by  fermio¬ 
nic  excitations  and  suggest  some  ways  the  latter  may  inform  the 
study  of  interacting  TIs  in  three  dimensions  [29].  We  close  with  a 
discussion  of  a  newly  emerging  direction  in  transition  metal  oxides: 
interaction-driven  topological  phases  at  interfaces  [30-33]. 

Our  paper  is  organized  as  follows.  In  Section  2  we  set  the 
convention  for  our  use  of  the  term  “topological  order”  in  this  paper 
and  provide  some  important  general  background,  including  Tl-QSL 
connections  in  one-dimension  and  important  results  for  interacting 
systems.  In  Section  3  we  detail  our  main  results  for  two-dimen¬ 
sions,  and  in  Section  4  we  describe  our  chief  three-dimensional 
results.  Einally,  in  Section  5  we  summarize  our  main  results  and 
give  our  view  of  interesting  future  directions  of  study. 


2.  Preliminaries 

2.3.  Topological  terminology 

Our  main  goal  in  this  paper  is  to  provide  examples  of  how  TIs 
and  QSLs  are  related  to  each  other,  and  how  the  study  of  each  can 
help  deepen  our  understanding  of  the  other.  An  important 
“connection”  between  these  two  types  of  phases  is  via  their 
topological  [34]  properties.  In  Section  3,  we  will  be  more  precise 
by  what  we  mean  by  “connection”.  We  note  here  that  it  does  not 
simply  mean  adiabatic  connection  of  the  electronic  states  because 
TIs  and  QSLs  are  distinct  phases  of  matter.  In  order  to  help  avoid 
confusion  later,  we  emphasize  that  we  will  use  “topological”  in 
the  broadest  sense  of  the  word:  there  is  some  mathematical 
invariant  characterizing  the  global  properties  of  the  system  that 
cannot  change  unless  a  bulk  gap  in  the  energy  spectrum  closes. 
Those  with  a  background  in  the  fractional  quantum  Hall  effect  or 
exotic  quantum  spin  liquids  commonly  take  “topological  order” 
to  be  synonymous  with  “non-trivial  ground  state  degeneracy” 
and  “non-trivial  (i.e.  fractional)  quantum  numbers/statistics”  of 
excitations  [5].  Since  a  large  part  of  this  paper  deals  with  TIs 
which  do  not  have  a  non-trivial  ground  state  degeneracy,  but 
nevertheless  do  have  topological  properties,  we  choose  a  more 
encompassing  use  of  the  terminology. 

Eor  purposes  of  this  paper,  we  would  say  that  Z2  time-reversal 
invariant  TIs,  integer  quantum  Hall  systems,  fractional  quantum 
Hall  systems,  and  QSLs  described  at  low  energies  by  a  Z2  gauge 
theory  all  possess  “topological  order”.  We  find  this  use  of 
terminology  convenient  for  highlighting  the  different  roles  that 
topological  invariants  and  conventional  order  parameters  play  in 
describing  a  state  of  matter.  This  way,  the  quantities  characteriz¬ 
ing  a  phase  are  divided  into  “non-local”  (topological)  and  “local” 
(conventional  order  parameter)  types.  (As  an  aside,  we  note  that 
such  divisions  can  be  subtle  and  sometimes  go  against  conven¬ 
tional  wisdom.  Eor  example,  superconductors  should  be  classified 
as  topologically  ordered  in  the  sense  of  Wen  [5],  as  they  have  a 
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non-trivial  ground  state  degeneracy  and  possess  no  gauge  invar¬ 
iant  local  order  parameter  [35].) 

2.2.  General  remarks  on  one-dimensional  systems 

In  this  topical  review,  we  are  mainly  concerned  with  two  and 
three-dimensional  systems.  However,  a  few  remarks  about  one¬ 
dimensional  systems  are  in  order.  One-dimensional  Fermi  sys¬ 
tems  are  special  for  a  number  of  reasons:  (i)  They  are  strongly 
correlated  at  low-energies  regardless  of  the  strength  of  the 
interactions,  and  they  possess  interesting  properties  such  as 
spin-charge  separation,  a  kind  of  fractionalization  of  fermion 
quantum  numbers  [36,37],  (ii)  Their  theoretical  study  is  amenable 
to  powerful,  non-perturbative  techniques  based  on  Abelian  and 
non-Abelian  bosonization  techniques  [36,37].  A  well-known  gap¬ 
less  phase  of  one-dimensional  fermions  conveniently  described 
with  these  techniques  is  the  Luttinger  liquid  [38].  In  recent  years 
it  has  become  clear  that  interesting,  highly  universal  “spin- 
incoherent”  Luttinger  liquid  regimes  also  appear  [39-41]. 

Focusing  our  attention  on  gapped  one-dimensional  systems,  a 
few  results  are  suggestive  of  possible  connections  between  TIs 
and  QSLs.  While  there  are  no  Z2  topological  insulators  in  one 
dimension  [8,42,43],  the  more  general  Altland  and  Zirnbauer 
classifications  [44,45]  of  non-interacting  topological  states  has 
provided  us  with  important  results  for  other  topological  phases 
resulting  from  different  sets  of  symmetries.  Of  particular  interest 
is  the  class  BDI  (time-reversal  symmetry  with  sublattice  symme¬ 
try  and  integer  spin),  which  has  a  Z  classification  in  one  dimen¬ 
sion  [42,43].  Fidkowski  and  Kitaev  [46]  have  recently  shown  in  a 
one-dimensional  model  in  class  BDI  that  the  Z  classification 
breaks  down  to  Zg  in  the  presence  of  interactions,  and  have 
further  generalized  these  results  to  other  symmetry  classes  [47]. 
In  other  words,  phases  that  are  topologically  distinct  at  the  non¬ 
interacting  level  can  be  adiabatically  connected  with  interactions 
that  are  slowly  turned  on  and  off  in  a  particular  way.  This 
example  proves  that  interactions  can  lead  to  important  modifica¬ 
tions  of  the  non-interacting  topological  classifications  [48].  A 
similar  breakdown  may  occur  in  two  and  three-dimensional 
systems,  although  to  our  knowledge  there  is  not  yet  a  concrete 
example. 

2.3.  Entanglement  as  a  tool  to  study  topological  order 

Because  topological  order  is  a  global,  non-local  property,  it  is 
rather  difficult  to  measure  in  the  general  situation.  (Most  physical 
responses  are  described  by  the  local  coupling  of  an  order  parameter 
to  an  external  perturbation,  or  generalized  “force”  [49].)  In  a  few 
important  cases,  such  as  the  integer  and  fractional  quantum  Hall 
effects,  and  the  Z2  topological  insulators,  the  topological  properties 
give  rise  to  certain  quantized  responses  [5,8,50,51].  These  examples 
have  topologically  protected  gapless  boundary  modes  that  dominate 
the  responses.  However,  often  there  is  no  obvious  “nice”  response 
that  can  be  computed  or  observed  in  experiment.  The  challenge  of 
understanding  the  many-body  quantum  mechanics  in  these  cases 
has  led  to  important  advances  in  the  understanding  of  quantum 
entanglement  in  recent  years  [52,53]. 

In  particular,  the  entanglement  entropy  [54-56]  and  the  entan¬ 
glement  spectrum  [57]  have  emerged  as  two  important  measures  of 
the  quantum  entanglement  and  the  topological  properties.  Typi¬ 
cally,  a  reduced  density  matrix  of  the  ground  state  is  computed  by 
tracing  over  some  degrees  of  freedom  (usually  real-space  coordi¬ 
nates  for  part  of  the  system).  The  topological  entanglement  entropy, 
computed  by  tracing  over  the  reduced  density  matrix  times  the 
logarithm  of  the  reduced  density  matrix,  contains  important  infor¬ 
mation  about  the  “quantum  dimensions”  of  the  excitations  in  the 
system  [55,56],  This  measure  is  useful  for  topological  systems  with 


excitations  that  have  “exotic”  quantum  numbers,  such  as  the 
fractional  quantum  Hall  effect  [58],  (The  non-interacting  22  topolo¬ 
gical  insulators  do  not  have  “exotic”  quantum  numbers,  so  the 
topological  entanglement  entropy  is  not  a  particularly  revealing 
quantity  in  this  case.  However,  in  a  recently  identified  strongly 
correlated  QSH*  phase  not  adiabatically  connected  to  the  QSH 
phase,  the  entanglement  entropy  can  serve  as  a  means  to  distin¬ 
guish  them  [59].) 

On  the  other  hand,  the  entanglement  spectrum  (ES)  [57]  has 
proven  to  be  useful  in  the  study  of  topological  insulators 
[22,60-64].  The  ES  is  found  by  writing  the  reduced  density  matrix 
as  an  exponential  of  a  fictitious  Hamiltonian.  The  eigenvalues  of 
this  fictitious  Hamiltonian  constitute  the  entanglement  spectrum. 
It  is  remarkable  that  the  ES,  which  is  determined  only  from  the 
ground-state  wave  function,  contains  important  information 
about  the  excitations  of  the  system  [57],  In  particular,  for  non¬ 
interacting  models  (such  as  Z2  topological  insulators  and  super¬ 
conductors)  it  has  been  shown  that  gapless  boundary  modes 
imply  a  gapless  ES  [60,61,65].  Moreover,  it  has  recently  been 
proved  by  Chandran  et  al.  [66]  that  in  fractional  quantum  Hall 
states  the  low-energy  modes  of  the  entanglement  spectrum  and 
edge-state  excitation  spectrum  are  in  one-to-one  correspondence. 
Qi  et  al.  have  obtained  similar  results  [67]. 

One  important  result  to  emerge  from  the  study  of  the  ES  of 
topological  insulators  is  that  a  gapless  ES  can  persist  under  some 
conditions  where  the  physical  edge  spectrum  becomes  gapped. 
For  example,  applying  a  magnetic  field  to  an  inversion  symmetric 
Z2  topological  insulator  will  gap  the  surface  states  but  leave  the 
entanglement  spectrum  gapless  [61,62].  The  same  is  true  if  time- 
reversal  symmetry  is  broken  in  more  subtle  ways  [62].  Strong 
interactions  that  gap  the  edge  modes  [68,69]  will  also  leave  a 
degeneracy  in  the  ES  [70],  Thus,  the  ES  “remembers”  the  under¬ 
lying  state  is  topological.  In  this  case,  it  turns  out  that  the  gapless 
nature  of  the  ES  is  protected  by  an  inversion  symmetry  [60-63], 
In  the  case  of  three-dimensional  insulators  with  inversion  sym¬ 
metry,  a  gapless  entanglement  spectrum  implies  an  “E  ■  B"  term 
in  the  action  even  if  the  boundary  contains  no  gapless  modes  [61  [.  In 
this  sense,  the  entanglement  spectrum  can  be  used  to  identify  a 
phase  of  matter  [63].  Recently  it  has  been  shown  that  the  break¬ 
down  of  the  Z  classification  (to  Zg)  of  class  BDI  in  the  presence  of 
interactions  [46,47]  can  be  understood  via  the  ES  [71], 

We  close  this  subsection  on  entanglement  by  returning  once 
more  to  one  dimension  and  QSLs.  Eor  aficianados  of  one-dimen¬ 
sional  spin  chains,  Haldane’s  classification  [72]  of  antiferromagne- 
tically  coupled  Heisenberg  spins  is  celebrated:  half-odd  integer 
multiple  spin  chains  are  gapless  (Luttinger  liquids)  while  integer 
spin  chains  are  gapped  and  possess  no  long-range  order,  i.e.  they  are 
gapped  QSLs.  Technically,  the  two  cases  are  distinguished  by  the 
Berry  phase  (topological)  term  in  the  effective  low-energy  action. 
Remarkably,  the  gapped  integer  spin  systems  contain  topological 
order  which  manifests  itself  as  a  non-local  string  order  parameter 
with  long  range  order  ]73],  and  protected  gapless  boundary  modes 
[74]  (similar  to  TIs).  In  fact,  ES  studies  of  integer  spin  chains  show 
phenomenology  quite  similar  to  TIs.  For  example,  the  physical  edge 
state  degeneracies  lead  to  the  same  degeneracies  in  the  ES  [75], 
Additionally,  perturbations  that  destroy  long-range  order  in  the 
string  parameter  and  edge-state  degeneracy  leave  the  ES  degeneracy 
in  tact,  so  long  as  inversion  symmetry  is  preserved  [75].  The  ES 
results  show  that  the  perturbations  to  spin  chains  that  destroy  the 
conventional  measures  of  the  topological  order  in  the  Haldane 
phase,  long-range  string  order  and  edge-state  degeneracies,  do  not 
destroy  the  Haldane  phase  itself  (so  long  as  the  bulk  gap  remains 
open).  (Even  in  gapless  one-dimensional  systems  the  ES  may 
provide  useful  information  [76].) 

Thus,  the  entanglement  properties  have  emerged  as  a  robust 
method  of  identifying  a  quantum  phase  [77],  It  is  remarkable  that 
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the  behaviors  of  the  ES  spectrum  for  TIs  and  QSLs  exhibit  such 
similar  phenomenology,  both  with  and  without  perturbations 
that  destroy  the  edge  modes.  As  we  will  see  in  the  next  section, 
there  are  more  direct  ways  that  TIs  and  QSLs  can  be  related  in 
higher  spatial  dimensions. 


3.  Two-dimensional  systems 

In  this  section  we  will  describe  a  class  of  non-interacting  and 
interacting  lattice  models  with  topological  order  in  two  spatial 
dimensions.  Some  of  these  are  spin  models  of  QSLs  and  some  are 
tight-binding  models  of  TIs.  We  will  show  that  under  a  certain  set  of 
conditions  there  is  a  precise  topological  connection  between  them. 

3.1.  Orientation 

The  first  prediction  of  a  quantum  spin  Hall  state  (a  two-dimen¬ 
sional  Tl)  in  a  specific  material  was  made  by  Bernevig  et  al.  [78]  for 
HgTe  quantum  wells.  This  prediction  was  soon  verified  experimen¬ 
tally  in  a  series  of  beautiful  experiments  in  Wurzburg  [79,80]. 

In  this  topical  review,  we  will  focus  on  models  inspired  by  the 
earlier  pioneering  work  of  Kane  and  Mele  [81,82],  which  first 
established  the  quantum  spin  Hall  state  as  a  new  topological 
phase  in  a  simple  tight-binding  model  on  the  honeycomb  lattice 
that  preserves  time-reversal  symmetry.  Their  Hamiltonian  is 

Wk-M  =  —  t  ^  (C?jCjo.-t-H.C.)  +  WcDW  +  WR 
(ijy.o- 

-hiAso  (*1) 

<<!/>>.“./* 

where  t  is  a  real  first-neighbor  hopping,  cj^  creates  a  particle  of 
spin  (7  on  site  i  and  annihilates  the  same  particle,  /tso  is  the 
(real)  parameter  representing  the  strength  of  the  second-neighbor 
spin-orbit  coupling,  e,j  =  (dj  x  dfj) / 1  dj  x  d?  |  is  a  vector  normal  to 
the  x-y  plane  describing  how  the  path  <  <i;>>  was  traversed, 
that  is,  either  “bending”  to  the  right  or  left  [81,82],  Here,  TYcdw 
and  Hr  are  charge  density  wave  and  Rashba  spin-orbit  terms, 
respectively.  The  CDW  term  describes  a  mean-field  staggered 
sublattice  potential  that  would  gap  the  Dirac  points  if  /tso  =  0  and 
put  the  model  in  a  trivial  insulator  phase.  The  Rashba  spin-orbit 
term  breaks  conservation  and  is  generally  present  if  the 
physical  system  has  broken  inversion  symmetry  with  respect  to 
the  plane  of  the  honeycomb  lattice. 

A  key  precedent  to  Kane  and  Mele’s  work  was  a  result  from 
Haldane  ]83]  in  which  it  was  shown  that  a  tight-binding  model  of 
spinless  fermions  on  the  honeycomb  lattice  with  broken  time-reversal 
symmetry  but  zero  net  magnetic  field  could  exhibit  a  quantum  Hall 
effect.  (In  other  words,  it  is  the  breaking  of  time-reversal  symmetiy 
that  is  key  for  the  quantum  Hall  effect,  not  the  net  magnetic  flux.) 
Haldane’s  model  is 

WH  =  -tX](c/Q  +  H.c.)-i-?YcDw  +  >t2  sgn(eij)cjcj,  (2) 

<!J>  <<&■>> 

where  t2  is  the  (real)  parameter  describing  the  string  of  the  second 
neighbor  hopping  and  sgn(eij)  is  + 1  depending  on  whether  the 
spinless  fermion  made  a  right  or  left  “turn”  along  the  path  <  <  y  >  > . 

As  Kane  and  Mele  emphasized  in  their  work  ]81,82],  in  a 
quantum  spin  Hall  system  that  preserves  quantization  \Hr  =  0 
in  Eq.  (1)],  the  state  can  be  represented  as  “two  copies”  of 
Haldane’s  model  (2) — one  for  spin  up  and  one  for  spin 
down — each  of  which  sees  an  opposite  sign  of  the  effective 
magnetic  field  and  therefore  preserves  time-reversal  symmetry 
overall.  The  presence  of  a  small  but  finite  Rashba  coupling  does 
not  change  the  phase;  it  only  destroys  the  conservation. 


While  the  topological  classification  of  the  quantum  Hall  effect  is 
based  on  the  Chern  number  [58]  (which  takes  the  integers  Z)  and 
the  two-dimensional  time-reversal  invariant  TIs  are  classified  by  Z2 
[1-3],  the  Haldane- Kane-Mele  correspondence  shows  that  the 
integer  quantum  Hall  effect  in  a  lattice  model  is  intimately  related 
to  the  quantum  spin  Hall  state.  In  particular,  if  a  spin  polarized 
tight-binding  model  realizes  a  quantum  Hall  effect,  combining  two 
time-reversed  copies  of  it  with  spin  degrees  of  freedom  is  guaran¬ 
teed  to  produce  a  quantum  spin  Hall  effect.  This  simple  observation 
is  actually  a  powerful  tool  for  generating  desired  topological  states  in 
two-dimensional  lattice  models.  We  will  use  this  approach  in 
several  different  ways  below,  including  generating  TIs  from  QSLs 
[62,84],  generating  QSLs  from  TIs  [18]  (tbe  Haldane-Kane-Mele 
correspondence  in  reverse),  and  generating  fractional  TIs  from  the 
fractional  quantum  Hall  effect  in  lattice  models  [85], 

3.2.  "Non-interacting"  lattice  models 

After  the  work  of  Kane  and  Mele  on  the  honeycomb  lattice 
[81,82]  there  was  a  flurry  of  activity  exploring  simple  non-interact¬ 
ing  tight-binding  lattice  models  on  a  variety  of  different  lattices  in 
two  and  three  spatial  dimensions.  For  example,  in  two  dimensions 
the  decorated  honeycomb  lattice  [84],  the  checkerboard  lattice  [86], 
the  square-octagon  lattice  [62],  the  kagome  lattice  ]87],  the  ruby 
lattice  [85],  and  others  [88]  have  been  shown  to  support  a  TI  phase 
in  simple  s-band  models.  If  multiple  orbitals  on  each  site  are  used, 
even  the  simple  square  lattice  can  support  a  TI  [78], 

A  natural  question  to  ask  is  “What  properties  of  a  lattice  model 
would  allow  it  to  support  a  TI  phase?”  One  way  to  view  this 
question  in  light  of  the  discussion  above  is  “What  properties  of  a 
spinless  fermion  lattice  model  would  allow  it  to  support  an  integer 
quantum  Hall  effect?”  When  cast  this  way,  there  are  immediate 
connections  to  so-called  “chiral”  QSLs,  a  class  of  spin  liquids  with 
broken  time-reversal  symmetry  and  a  thermal  Hall  conductance.  In 
fact,  one  of  the  early  models  relating  spin  liquids  to  the  (fractional) 
quantum  Hall  effect  was  proposed  by  Laughlin  himself  [89], 

In  recent  years,  an  important  class  of  exactly  solvable  spin 
models  known  as  “Kitaev”  models  has  emerged  ]10,n].  These 
models  can  be  used  to  generate  concrete  examples  of  various 
types  of  QSLs  (previous  studies  often  relied  on  uncontrolled 
approximations).  In  particular,  chiral  spin  liquids  with  a  quan¬ 
tized  thermal  Hall  conductance  can  be  found  [18,90,91].  As  we 
now  show,  such  models  can  be  adiabatically  deformed  to  an 
integer  quantum  Hall  “Haldane-type”  model,  which  in  turn  can  be 
used  to  construct  a  two-dimensional  TI  via  the  Haldane-Kane- 
Mele  correspondence  described  above. 

As  a  simple  example,  consider  the  Kitaev  model  on  the 
honeycomb  lattice  ]11], 

WlCitaev  =  (3) 

Oi),  <ij>y  <V>z 

where  the  spin-1 /2  components  couple  to  each  other 

along  the  x,y,z-labeled  links,  respectively.  The  symmetry  of  the 
Hamiltonian  (3)  is  unusual  (not  Heisenberg,  XY,  or  Ising,  for 
example),  but  there  are  recent  proposals  suggesting  it  would 
arise  naturally  in  some  solid  state  contexts  [92,93],  The  main 
reason  for  choosing  the  form  (3)  is  that  it  admits  an  exact  solution 
with  a  spin-liquid  ground  state  [11]. 

The  Hamiltonian  (3)  is  easily  solved  by  choosing  an  alternate 
representation  of  the  spin  degrees  of  freedom:  One  can  write  the 
spins  as 

af  =  itijt]f,  a^.=iif.tfl,  =  (4) 

where  the  11° ,ri’j  ,tfj  .ifj  are  four  “flavors”  of  Majorana  fermions  on 
the  site  j.  They  have  the  properties  =  ((7“)^  =  1  on  each 

site  for  a  =  0,x,y,z,  and  they  satisfy  the  anticommutation  relations 
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=  ISabSij.  They  contain  “half  the  degrees  of  freedom  of  a 
“normal"  fermion  (one  that  satisfies  the  usual  anticommutation 
relations).  Importantly,  the  transformation  (4)  preserves  the  spin 
commutation  relations  [o'^.oj’]  =  2i£“*'‘^oT  provided  we  impose  the 
constraint  Dj  =  =  1  on  each  site.  (The  Majorana  repre¬ 

sentation  (4)  actually  doubles  the  Hilbert  space,  so  one  needs  to 
project  to  the  physical  Hilbert  space.)  Under  the  representation 
(4)  the  Kitaev  Hamiltonian  is  transformed  to 

Wrataev  =  (^> 

oj>x  (ijy,  <v>z 

where  Ujj  =  -irjfrjJ  is  a  Z2  gauge  field  with  eigenvalues  + 1 
defined  along  an  a-link.  The  product  of  six  different  u,j  around 
hexagonal  plaquettes  on  the  honeycomb  lattice  defines  a  Z2  gauge 
flux  through  that  plaquette,  Wp  =  H  Uij-  usefulness  of  the 

representation  (4)  is  that  now 

[Uy,Ufc/]  =  [Uy.TfKitaev]  =  [klfp»^Kitaev]  =  (6) 

SO  that  the  Uy  are  non-dynamical  gauge  fields  and  Wp  are  non- 
dynamical  gauge  fluxes.  In  other  words,  the  Uy  in  Eq.  (5)  can  be 
taken  to  be  constants,  rather  than  operators  where  the  value  of 
the  constants  are  determined  (up  to  gauge  equivalent  choices)  by 
the  fluxes  Wp  around  the  plaquettes  in  the  lattice.  For  a  bipartite 
lattice  (like  the  honeycomb  lattice)  it  is  known  that  the  uniform 
flux  configuration  is  lowest  in  energy  [94],  and  for  non-bipartite 
lattices  detailed  numerical  studies  indicate  the  same  [18,90]. 

With  the  Ujj  as  constants,  it  is  clear  that  the  interacting  spin 
Hamiltonian  (3)  has  now  been  reduced  to  a  non-interacting 
hopping  Hamiltonian  (5).  (Recall  that  =  Of  course,  a 

non-interacting  Hamiltonian  can  always  be  solved  exactly.  In 
the  present  case,  it  is  useful  to  Fourier  transform  to  momentum 
space  as  one  would  for  any  tight-binding  Hamiltonian  and 
determine  the  energy  bands  E(k).  Because  of  the  self-adjoint 
property  of  the  Majorana  fermions,  half  of  the  states  are  redun¬ 
dant  and  E(k)  =  -E(-k).  One  “unusual”  property  of  the  tight- 
binding  model  (5)  is  that  the  hopping  amplitudes,  ijuy,  are  purely 
imaginary  (assuming  y  is  real). 

The  procedure  for  solving  any  Kiteav-type  model  on  any  lattice 
is  exactly  the  same  as  above,  even  if  the  lattice  is  different  and/or 
the  spin  is  larger  [18,95-99].  Staying  with  the  honeycomb  lattice 
model  for  the  moment,  it  is  worth  emphasizing  that  Kitaev’s 
solution  also  determines  whether  the  ground  state  of  the  spin 
system  spontaneously  breaks  time-reversal  symmetry  [11].  The 
statement  is  the  following:  If  the  lattice  contains  plaquettes  with 
only  an  even  number  of  sides  (such  as  the  honeycomb  lattice)  the 
ground  state  will  not  break  time-reversal  symmetry.  If  there  are 
any  plaquettes  with  an  odd  number  of  links  (such  as  the 
decorated  honeycomb  lattice  [90]  or  kagome  lattice  [18])  the 
ground  state  will  spontaneously  break  time-reversal  symmetry 
[11[.  This  statement  follows  because  time-reversal  symmetry 
sends  Uy^-Uy,  which  will  change  the  sign  of  Wp  on  odd-sided 
plaquettes,  but  not  even-sided  ones. 

Applying  this  to  the  honeycomb  lattice  model  (5)  we  can 
immediately  infer  that  since  the  model  does  not  break  time- 
reversal  symmetry,  the  band  structure  must  have  zero  Chern 
number.  On  the  other  hand,  the  addition  of  the  term  [12] 

H'=  ^ (7) 

(ijk)  e  A  (ijk)  e  V 

explicitly  breaks  time-reversal  symmetry  and  leads  to  a  second- 
neighbor  hopping  term 

W'=  -if^inkh°tl°k+  VuyUjf,ri°h°k-  (8) 

<<i(t>>sA  (<ik))eV 

where  A  and  V  denote  up  and  down  pointing  triangles  of  (ijk), 
with  i  always  taken  to  be  the  left  most  point.  Remarkably,  if  the 


Fig.  1.  (Color  online)  The  decorated  honeycomb  (sometimes  called  the  “star") 
lattice.  An  exactly  solvable  chiral  spin  liquid  model  with  a  finite  Chern  number 
and  non-Abelian  anions  has  been  found  on  this  lattice  [90].  By  the  Kitaev- 
Haldane-Kane-Mele  correspondence,  this  implies  a  topological  insulator  phase  in 
a  simple  s-band  model  also  exists  on  this  lattice  [84]. 

Majorana  band  structure  of  the  Hamiltonian  (5)  with  Eq.  (8)  are 
taken  to  represent  spinless  fermions  (as  opposed  to  Majorana 
fermions),  the  states  can  be  adiabatically  deformed  into  the 
Haldane  model  (2)  (which  contains  purely  real  first-neighbor 
hopping  and  imaginary  second-neighbor  hopping)  on  the  honey¬ 
comb  lattice  [12[.  This  implies  that  the  Hamiltonian  (5)  with 
Eq.  (8)  is  a  topological  chiral  spin  liquid  with  Chern  number  +  1 ! 

This  result  establishes  a  Kitaev-Haldane-Kane-Mele  correspon¬ 
dence  on  the  honeycomb  lattice  that  generalizes  the  earlier  Haldane- 
Kane-Mele  correspondence  to  include  chiral  spin  liquids  with  finite 
Chern  numbers  (provided  the  chiral  spin  liquids  can  be  described  by 
fermions  with  a  band  structure — there  is  not  such  an  obvious 
connected  to  spin  liquids  described  by  a  bosonic  representation  of 
the  spin  degrees  of  freedom).  As  we  now  show,  the  presence  of  a 
chiral  spin  liquid  with  finite  Chern  number  on  the  honeycomb  lattice 
[12]  when  there  is  also  an  integer  quantum  Hall  effect  [83]  and  a  TI 
phase  [81,82]  is  a  special  case  of  a  more  general  result. 

The  first  exactly  solvable  model  of  a  chiral  spin  liquid  with  a 
finite  Chern  number  that  included  only  two-spin  interactions 
was  proposed  by  Yao  and  Kivelson  [90]  on  the  decorated 
honeycomb  lattice  shown  in  Fig.  1.  The  Hamiltonian  is  of  the 
Kitaev  type  and  is  given  by 

<.ijyx  <(/>,  <(/>z 

+  (9) 

<y>v  Oi>y  <ij>2' 

where  the  exchange  is  given  by  J  along  the  links  of  the  triangles 
and  /  along  the  links  between  triangles,  as  shown  in  Fig.  1.  The 
Hamiltonian  (9)  is  solved  using  the  transformation  (4)  and  a 
fermion  hopping  model  with  purely  imaginary  hopping  para¬ 
meter  results,  similar  to  what  is  found  in  Eq.  (5).  Because  of  the 
triangular  plaquettes  on  the  lattice,  the  ground  state  sponta¬ 
neously  breaks  time-reversal  symmetry.  For  0  <  j'/j  <  V3,  the 
ground  state  is  a  gapped  chiral  spin  liquid  with  Chern  number 
+  1,  while  for  J'/J  >  \/3,  the  ground  state  is  gapped  chiral  spin 
liquid  with  zero  Chern  number  [90].  If  the  Kitaev-Haldane- 
Kane-Mele  correspondence  described  above  on  the  honeycomb 
lattice  is  to  prove  general,  Yao  and  Kivelson’s  discovery  of  the 
chiral  spin  liquid  with  finite  Chern  number  in  a  Kitaev  model 
should  imply  that  there  is  a  Tl  in  an  s-band  model  on  the 
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Fig.  2.  The  decorated  honeycomb  lattice  energy  bands  along  high  symmetry 
directions  in  the  Brillouin  zone  (from  Ref.  [84]).  In  (a)  Iso  =  0  and  in  (b)  Aso  =  0.1 1. 
There  are  Dirac  points  at  K  and  K'  (not  shown)  and  quadratic  band  crossing  points 
(QBCP)  at  r  in  (a),  while  in  (b)  Aso  ¥=  0  opens  up  a  gap  at  each  of  these  points  and 
destroys  the  flat  bands. 
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Fig.  3.  Phase  diagrams  for  the  decorated  honeycomb  lattice  with  t  and  t'  real  in 
the  absence  of  a  staggered  on-site  potential  and  no  Rashba  coupling  (from  Ref. 
[84]).  Several  filling  fractions / are  shown  (lower  left  corner).  For  fixed/and  Aso  it 
is  possible  to  drive  a  transition  between  a  topological  insulator  and  a  non- 
topological  phase  by  varying  the  ratio  t'/t. 


that  at  all  filling  fractions  there  is  a  finite  region  of  parameter  space 
occupied  by  a  T1  (QSH)  phase. 

We  now  show  that  the  T1  phase  at  1  /2  filling  is  “connected”  to 
the  chiral  spin  liquid  phase  of  Eq.  (9)  via  the  Kitaev-Haldane- 
Kane-Mele  correspondence  demonstrated  on  the  honeycomb  lat¬ 
tice  by  Lee  et  al.  [12].  We  do  this  by  showing  that  the  Hamiltonian 
describing  the  band  structure  of  Eq.  (9)  (with  “normal”  spinless 
fermions  replacing  the  Majorana  fermions)  with  only  purely 
imaginary  first-neighbor  hopping  can  be  adiabatically  deformed 
to  one  spin  component  of  Eq.  (10)  (with  Hcdvj  =  'Hr  =  0)  which  has 
purely  real  first-neighbor  hopping  and  purely  imaginary  second- 
neighbor  hopping.  This  is  done  by  applying  a  combination  of  gauge 
transformations  and  adiabatic  tuning  of  the  real  and  imaginary 
components  of  first  and  second  neighbor  hopping  parameters.  The 


Fig.  4.  (Color  online)  (from  Ref  [84])  Schematic  illustration  of  the  continuous  path 
which  adiabatically  connects  the  representative  free  fermion  model  of  the  Kitaev 
model  on  the  decorated  Honeycomb  lattice  with  the  spinless  model  with  real  t,  t' 
and  Aso.  The  patterns  (o)  and  (1)  are  equivalent  and  correspond  to  Kitaev's  model 
whereas  (v)  represents  the  model  with  real  t,  t'  and  Asq.  The  adiabatic  deformation 
does  not  lead  to  a  gap  closing  and  the  Chern  number  stays  constant.  See  Fig.  5.  This 
establishes  the  topological  connection  between  the  two  models. 


decorated  honeycomb  lattice.  We  have  shown  this  is  indeed  the 
case  [84]. 

Consider  the  s-band  Hamiltonian  on  the  decorated  honeycomb 
lattice, 

Wdhc-TB  =  -t  ^  clcj„-t'  + 

+WcDw+WR  +  i2.so  ^  ey  -Scipcj^Cj/i,  (10) 

<  (i]}  >.a,/J 

which  describes  nearest-neighbor  hopping  on  the  triangles  “A” 
with  amplitude  t  and  between  triangles  “A  -►  A”  with  amplitude  t', 
and  the  remaining  terms  describe  the  same  physics  as  represented 
in  Eq.  (1 ).  The  bands  for  "Hcdw  =  Wr  =  0  are  shown  in  Fig.  2,  and  the 
phase  diagrams  at  various  filling  fractions  are  shown  in  Fig.  3.  Note 


-1  -0.5  0  0.5  1 

</>/t 


Fig.  5.  (Color  online)  Contour  plot  of  (a)  the  Chern  number  and  (b)  the  gap  of  the 
model  as  a  function  of  (^,  Aso )  at  half-filling  for  J  =  f.  Here  (j)  is  the  phase  angle  on  the 
nearest  neighbor  hopping  parameters  given  in  Ref  [84].  Also  shown  is  a  possible  path 
which  adiabatically  connects  the  flux  patterns  (iii)-(v)  defined  in  Fig.  4. 


G.A  Fiete  et  al.  /  Physica  E  44  (2012)  845-859 


851 


procedure  is  shown  in  Fig.  4,  and  the  panels  in  Fig.  5  explicitly 
show  that  the  Chern  number  remains  unchanged  and  the  gap 
remains  open  in  this  process,  establishing  the  topological  equiva¬ 
lence,  and  the  Kitaev-Haldane-Kane-Mele  correspondence  on  the 
decorated  honeycomb  lattice  [84], 

It  should  be  clear  from  the  examples  of  the  decorated 
honeycomb  lattice  and  the  honeycomb  lattice  above  that  the 
Kitaev-Haldane-Kane-Mele  correspondence  is  rather  general. 
However,  it  is  also  true  that  one  can  invert  the  argument  to  use 
the  existence  of  a  TI  in  an  s-band  model  on  a  particular  lattice 
to  infer  a  chiral  spin  liquid  with  a  finite  Chern  number  in  a 
Kitaev-type  model.  Indeed,  we  have  used  the  work  of  Guo  and 
Franz  that  established  a  TI  on  the  kagome  lattice  [87]  to  motivate 
the  study  of  a  Kitaev-type  model  on  the  same  lattice  [18],  Not 
only  did  we  find  a  new  exactly  solvable  chiral  spin  liquid  model 
with  finite  Chern  number,  but  also  a  phase  that  possess  a  stable 
spin  Fermi  surface  [18],  We  note  that  unstable  spin  Fermi  surfaces 
have  been  reported  in  other  models  [95,100,101]. 

Because  the  kagome  lattice  has  sites  with  four  links  (as 
opposed  to  three  on  the  honeycomb,  decorated  honeycomb,  and 
decorated  square  lattice),  the  effective  spin  in  a  Kitaev-type 
model  must  be  at  least  3/2  [96].  These  higher-spin  Kitaev  models 
require  that  the  spin-1 /2  Pauli  spin  matrices  be  generalized  to 
Gamma  matrices  [95].  The  spin  Hamiltonian  we  study  [18]  on  the 
kagome  lattice  is 

^=Ja  E  rjrf+j^  E  rfrf+j,J2rf 

<lJ>eA  <ij>EV  i 

+h'  E  rrrf+h'  E  (11) 

<IJ>eA  <lJ>sV 

where  we  have  distinguished  the  nearest  neighbor  couplings  Jy  as 
Ja,  Jv  3nd  J'y  as  and  jy'  based  on  whether  the  link  <  ij)>  belongs 
to  an  up  (A)  triangle  or  down  (V)  triangle  (the  kagome  lattice  can 
be  viewed  as  a  corner  sharing  network  of  up  and  down  triangles), 
and  <  ij  >  is  taken  in  the  counterclockwise  sense  for  each  triangle. 
Locally,  the  five  T-matrices  satisfy  a  Clifford  algebra, 

{r[',rf)  =  2b“^  where  a,fa  =  l _ 5,  and  r""  ^  [r“,r'']/(2i).  In 

terms  of  the  components  of  the  spin  S=  3/2  operators  [95,102], 

=  r==(S")2-|.  (12) 

With  the  identification  (12),  it  is  clear  the  model  (11)  has  a  global 
Ising  spin  symmetry  under  180°  rotations  about  the  z-axis,  and 
possesses  time-reversal  symmetry  (TRS)  (although  TRS  will  be 
spontaneously  broken  in  the  ground  state  from  the  triangular 
plaquettes  as  we  described  earlier)  in  addition  to  the  translational 
and  threefold  rotational  lattice  symmetry  mentioned  above. 

In  order  to  solve  the  higher-spin  G-matrix  model  (11),  a 
Majorana  representation  similar  to  that  for  spin-1 /2,  (4),  is 
employed 

n  =  i^Ci.  rf  =  iCidi,  =  (13) 

with  a  =  1,2,3, 4.  There  are  thus  six  Majorana  species  on  each  site 
i:  The  Majorana  representation  enlarges  the 

spin-3/2  Hilbert  space,  so  that  one  must  enforce  the  constraint 
Dj  =  -r/rfrfrfrf  =  =  l,  namely  for  any  physical 

state  I W  >  phys  Di\Wy  p^yj  =  | '/'  >  phys  for  any  i.  The  remainder  of  the 
steps  are  similar  to  the  spin-i/2  version  and  one  obtains  a 
Majorana  tight-binding  model  on  the  kagome  with  imaginary 
hopping  parameters  [18].  Solving  the  model  on  a  strip  geometry 
in  the  gapped  phase  with  a  finite  Chern  number  (  +  2)  yields  the 
results  shown  in  Fig.  6.  The  gauge  dependent  but  stable  Fermi 
surface  in  a  gapless  phase  is  shown  in  Fig.  7. 


Fig.  6.  (a)  Band  structure  on  a  cylindrical  geometry  for  =J^'  =  '[  .0, 
jy  =  0.8,  Js  =  0.  There  are  two  gapless  chiral  Majorana  edge  states  (dotted 
lines)  which  overlap  on  each  other,  (b)  For  =  = 

and  75=0.1,  the  two  gapless  edge  states  separate.  These  ground  states  are  thus 
CSLs  with  a  spectrum  Chern  number  (  ±  2)  and  the  vortices  obey  Abelian  statistics. 
From  Ref.  [18]. 
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Fig.  7.  The  Fermi  surface  (solid  line)  for  the  flux  configuration  {7i/2,n/2,n}  and 
UaJvJaMvVsI  ^  {1.0, 0.3, 0.8, 0.5, 1.4}.  (See  Ref.  [18]  for  a  description  of  notation.) 
The  dashed  hexagon  is  the  Brillouin  zone  boundary.  Note  that  there  is  only  one 
Fermi  pocket  for  this  set  of  parameters  and  the  three  pockets  shown  are  related  by 
a  reciprocal  lattice  vector  and  are  thus  equivalent. 

The  results  of  this  subsection  have  demonstrated  the  Kitaev- 
Haldane-Kane-Mele  correspondence  and  illustrated  how  TIs  can  be 
generated  from  a  certain  class  of  QSLs,  and  vice-a-versa.  Thus,  we 
have  established  a  method  of  generating  certain  types  of  desired 
states.  Moreover,  along  the  way  we  have  also  discovered  new  QSLs, 
such  as  the  gapless  varieties  with  a  stable  spin  fermi  surface  [18].  In 
all  these  "non-interacting”  models,  the  topological  relationship 
between  the  TI  and  QSL  models  is  based  on  the  Chern  number  and 
an  underlying  band  structure  of  free  fermions.  In  the  next  subsection 
we  turn  to  the  study  of  interacting  models  in  two  dimensions. 
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As  a  final  remark,  we  note  there  are  also  deep  connections 
between  Kitaev  models,  the  quantum  Hall  effect,  topological 
insulators,  and  topological  superconductivity  [11,58,65,103-107], 
The  topological  Kitaev  models  [19]  and  the  TI  insulators  (including 
the  integer  quantum  Hall  effect)  [20-24]  show  that  disorder  can 
play  a  useful  role  in  stabilizing  topologically  ordered  states  and 
could  play  a  role  in  guiding  the  experimental  discovery  of  topolo¬ 
gical  quantum  spin  liquids  [19[.  The  non-interacting  models 
discussed  in  this  section  can  be  viewed  as  mean-field  states  in  a 
more  general  interacting  model. 

3.3.  Interacting  lattice  models 


ruby  [85[.  As  we  show  below,  the  cases  of  quadratic  band 
touching  points  often  have  topological  phases  as  the  leading 
instabilities  with  respect  to  interactions  [124].  This  effectively 
enlarges  the  candidate  systems  for  Tl  physics  to  include  those 
without  strong  spin-orbit  coupling  but  with  “fiat”  band  touching 
points  in  the  non-interacting  band  structure.  Another  notable 
example  is  bilayer  and  multilayer  graphene  [130]. 

The  kagome  lattice  serves  as  an  important  example  case.  Its 
non-interacting  band  structure  in  an  s-band  nearest  neighbor 
hopping  model  possess  both  Dirac  points  (1/3  filling)  and  quad¬ 
ratic  band  touching  points  (2/3  filling)  [87[.  We  study  the 
following  Hamiltonian  [124]: 


While  the  Kitaev  models  are  interacting  in  terms  of  the 
underlying  electronic  degrees  of  freedom  (local  moments),  we 
discussed  them  under  the  heading  of  “non-interacting”  lattice 
models  because  they  admit  a  free  (Majorana)  fermion  representa¬ 
tion  and  therefore  can  be  understood  in  terms  of  the  properties  of 
their  band  structure.  There  are  few  exactly  solvable  interacting 
models  in  two  dimensions,  and  those  that  exist  often  possess 
special  symmetries  [108].  Thus,  one  must  look  for  approximate 
methods  to  address  the  physics  of  interactions. 

For  a  certain  class  of  Hamiltonians  (those  without  a  fermion 
sign  problem),  quantum  Monte  Carlo  is  a  powerful  theoretical 
method.  Given  the  central  role  the  honeycomb  lattice  has  played 
in  the  understanding  of  Tls  and  QSLs,  it  is  natural  to  begin  an 
investigation  of  interacting  models  on  this  lattice.  In  the  last  year 
the  simplest  interacting  model,  the  Hubbard  model  (no  spin-orbit 
coupling  terms),  has  been  shown  to  possess  a  (apparently  gapped) 
spin-liquid  phase  at  half-filling  and  intermediate  coupling  [109]. 
Moreover,  this  spin  liquid  has  been  shown  to  be  robust  in  the 
presence  of  a  second-neighbor  spin-orbit  coupling  in  quantum 
Monte  Carlo  [110]  and  dynamical  mean-field  [111]  studies  of  the 
so-called  Kane-Mele-Hubbard  model.  Related  quantum  Monte 
Carlo  studies  on  bulk  and  edge  instabilities  in  the  Kane-Mele- 
Hubbard  model  have  been  reported  [111-114],  as  well  as  closely 
related  exact  diagonalization  studies  in  spinless  systems  with 
nearest-neighbor  interactions  for  small  system  sizes  [115]. 

Aside  from  bosonization  studies  of  interactions  in  the  one¬ 
dimensional  helical  liquids  that  exist  along  the  boundary  of  a  two- 
dimensional  Tl  [116-123],  most  studies  of  interactions  in  bulk 
have  relied  on  some  type  of  mean-field  theory  [17,59,86,124-127]. 
One  of  the  most  important  results  to  emerge  from  these  studies  is 
that  topological  insulators  can  be  spontaneously  generated  from 
interactions  at  the  mean-field  level  even  when  there  is  no  micro¬ 
scopic  spin-orbit  coupling  present.  This  was  first  explicitly  shown  by 
Raghu  et  al.  on  the  honeycomb  lattice  [17[.  In  model  without  any 
intrinsic  spin-orbit  coupling,  TI  phases  are  produced  by  the 
spontaneous  generation  of  spin-orbit  coupling  from  first  and 
second-neighbor  interactions  treated  at  the  mean-field  level 
[17,30,86,124-126]. 

The  stability  of  a  non-interacting  band  structure  to  interac¬ 
tions  depends  on  a  large  extent  on  the  form  of  the  low-energy 
dispersions  and  band  crossing  points.  For  example,  it  is  well- 
known  that  at  1/2  filling  on  the  honeycomb  lattice  the  non¬ 
interacting  dispersions  are  Dirac-like  [128].  Dirac  points  are 
perturbatively  stable  to  interactions  [17,86,88,124].  In  practice, 
this  implies  that  a  critical  interaction  strength  is  required  to 
drive  a  phase  transition  away  from  the  gapless  state  [86[. 
This  is  qualitatively  consistent  with  both  mean-field  studies 
[17,124,129]  and  quantum  Monte  Carlo  studies  [109].  By  contrast, 
two  bands  that  touch  quadratically  or  flatter  are  infinitesimally 
unstable  to  interactions  in  two  dimensions  [86,124].  Quadratic 
band  touching  points  are  known  to  occur  on  many  of  the  lattices 
that  support  Tls  in  simple  s-band  models:  decorated  honeycomb 
[84],  kagome  [87],  checkerboard  [86],  square-octagon  [62],  and 


^kagome  —  1  ^  ^  Cjfj  -1-  'Uj |  H,- 1 

<i.j)  i 

+  V]  ^  n,nj  +  V2^  riiUj  +  ^3^  mj  (14) 

<W>  «ij»  €'J'» 

at  the  mean-field  level.  Here,  c!p  annihilates  (creates)  a  fermion 
on  site  i  with  spin  cr=  I,!,  nirr  =  cJ^Cj^  and  ni  =  ^^nirj.  The  sums 
run  over  nearest-neighbor  <  ij  > ,  second-neighbor  <  i,j  > ,  or  third- 
neighbor  bonds  The  hopping  amplitude  is  denoted  by  t 

and  the  parameters  V^,  V2,  and  V3  quantify  the  nearest-neighbor, 
second-neighbor  and  third-neighbor  repulsion,  respectively. 

To  perform  the  mean-field  calculation,  we  decouple  the  on-site 
interaction  according  to  [124] 
nq  riii  ~  n,,  <  n,!  >  +  <  n,|  >  n,^  -  <  r!,|  >  <  114  > 

Cii  <  cj^  C4  >  -  <  cj^ C4  >  C4  -t  <  ct^  C4  >  <  C4  > .  (15) 


Fig.  8.  Filling  fraction  1/3  on  the  kagome  lattice  with  Dirac  points  in  the  non¬ 
interacting  band  structure.  The  U-V2  phase  diagram  for  Vi  =  0  and  V3  =  0.4t  in  Eq. 
(14).  An  interaction-driven  Tl  appears  for  finite  U  and  V2  and  requires  a  critical 
interaction  strength  to  appear,  and  some  fine-tuning  of  parameters.  Solid  lines 
indicate  first  order  and  dashed  lines  second  order  transitions.  From  Ref.  [124]. 


Fig.  9.  The  U-V]  phase  diagram  for  V2  =  V3  =  0  at  2/3  filling  fraction  (quadratic 
band  touching  point)  in  Eq.  (14).  Note  that  for  small  generic  interactions  the 
leading  instability  is  a  topological  phase,  either  TI  or  the  quantum  anomalous  Hall 
phase  (QAH).  From  Ref.  [124]. 
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Fig.  10.  (Color  online)  Schematic  of  a  drag  measurement  between  two  QSH 
systems  from  Ref.  [121].  A  current  Ii  is  driven  along  the  upper  edge  of  the  lower 
QSH  system  and  through  electron-electron  interactions  a  voltage  V2  is  induced  in 
the  lower  edge  of  the  upper  QSH  system.  A  magnetic  field  h  is  applied  in  the  plane 
of  wires,  perpendicular  to  the  spin  quantization  axis  (assumed  perpendicular  to 
the  plane  of  QSH  systems).  Time-reversed  Kramer’s  pairs  are  indicated  for  the  two 
edges.  A  QSH  on  top  of  QSH  geometry  could  also  be  used. 

We  assume  the  mean-field  solutions  are  described  by  a  collinear 
spin  alignment  and  therefore,  without  loss  of  generality,  we  set 
>  =  0  [124].  The  further-neighbor  interaction  is 
decoupled  in  a  similar  way  [124]: 

riirij  ^  Tii  <  Hj  >  +  <  n,  >  n,- <  n;  >  <  nj  >  -^  (clcjp  <  c1^c,„  > 

Ct/i 

+  <  cLovi  >  > )  ■  ('•6^ 

The  results  of  the  mean-field  calculations  are  shown  in  Figs.  8 
and  9.  For  a  detailed  description  of  the  phases  appearing  in  the 
phase  diagrams,  see  Ref.  [124]. 

The  important  message  we  would  like  to  emphasize  is  that  for 
Dirac  points,  Fig.  8,  topological  phases  require  a  critical  interac¬ 
tion  strength  that  can  be  rather  large,  and  appear  to  require  some 
amount  of  “fine  tuning”  [125],  For  example,  note  that  V]  =0  in 
Fig.  8.  On  the  other  hand,  when  quadratic  band  touching  points 
are  present  in  the  non-interacting  band  structure.  Fig.  9,  topolo¬ 
gical  phases  can  appear  as  the  leading  instabilities,  without  any 
“fine  tuning”  of  the  interactions.  Thus,  two-dimensional  band 
structures  with  flat  band  touching  points  but  small  (or  even 
vanishing)  spin-orbit  coupling  are  excellent  candidates  for  topo¬ 
logical  phases  [130].  General  arguments  and  functional  renorma¬ 
lization  group  studies  (that  treat  fluctuations)  suggest  the  mean- 
field  results  may  be  reliable  and  that  the  T1  phase  is  likely 
preferred  over  the  quantum  anomalous  Hall  (QAH)  state  when 
they  are  degenerate  at  the  mean-field  level  [17[. 

A  final  point  on  which  we  would  like  to  remark  concerning 
interacting  topological  insulators  in  two  dimensions  is  how 
interactions  may  help  identify  the  phase  itself.  For  example,  if 
two  edges  of  different  quantum  spin  Hall  systems  are  brought 
close  to  one  another  the  edges  will  interact  via  the  Coulomb 
interactions  of  the  electrons  in  each  edge  (See  Fig.  10).  Under 
some  circumstances,  this  can  lead  to  novel  quasi-one  dimensional 
phases  [122],  However,  if  the  interactions  are  not  too  strong  (the 
edges  are  far  enough  apart),  a  current  driven  in  one  edge  will 
result  in  a  voltage  drop  along  the  other  edge  in  the  so-called 
“Coulomb  drag”  effect  [131-137].  The  properties  of  the  edge  are 
reflected  in  the  temperature  and  magnetic  field  dependence  of 
the  drag.  For  identical  edges  described  by  a  helical  liquid  [68,69], 
the  Coulomb  drag  has  a  novel  magnetic  field  dependence  that  is 
“turned  on”  in  the  presence  of  a  time-reversal  symmetry  breaking 
magnetic  field  [121], 

(17) 

where  h  is  the  strength  of  the  applied  magnetic  field  and  K  is  the 
Luttinger  parameter  of  the  helical  edge  system  [121].  The  role  of 
the  finite  magnetic  field  is  to  open  a  back-scattering  channel  for 
electrons.  There  is  also  a  weak  field  dependence  of  the  Luttinger 
parameter  K  [121]. 


3.4.  Fractional  topological  insulators  and  flat  band  fractional 
quantum  Flail  effect 

The  discussion  in  Section  3.2  which  related  the  integer 
quantum  Hall  effect  to  TIs  immediately  raises  the  question  of 
whether  there  is  a  “fractional”  analog  to  TIs  that  depends  on  an 
intrinsic  way  on  electron-electron  interactions  in  the  same  way 
the  fractional  quantum  Hall  effect  does.  In  particular,  is  there  a 
time-reversal  invariant  topological  insulator  that  is  not  adiabati- 
cally  connected  to  a  free  electron  model?  Theoretically,  the 
answer  is  affirmative  [13-15]. 

In  two  dimensions,  we  are  guided  in  a  study  of  fractional 
topological  insulators  by  the  fractional  quantum  Hall  effect.  If  one 
has  in  hand  a  spinless  fermion  model  that  realizes  a  fractional 
quantum  Hall  effect,  an  conserving,  time-reversal  invariant 
model  can  be  straightforwardly  constructed  by  generalizing  the 
Haldane-Kane-Mele  correspondence  to  the  interacting  case 
[13,14].  Only  very  recently  has  the  explicit  construction  been 
understood  for  the  non-S^  conserving  case  (in  the  presence  of 
Rashba  spin-orbit  coupling,  for  example)  [15[. 

A  key  to  understanding  lattice  models  of  fractional  topological 
insulators  in  two  dimensions  is  thus  understanding  the  fractional 
quantum  Hall  effect  in  lattice  models.  (While  this  appears 
necessary,  it  is  not  sufficient  for  constructing  a  general  many- 
body  wavefunction  [15,138-141]).  For  a  while,  this  understand¬ 
ing  eluded  the  community  whose  picture  was  based  on  trial 
wavefunctions  and  effective  low-energy  theories  [58],  but  in  the 
past  year  important  progress  has  been  made. 

Attention  has  focused  on  a  class  of  insulators  with  nearly  flat 
bands  that  possess  finite  Chern  numbers.  With  a  finite  Chern 
number,  a  partial  filling  of  the  flat  bands  can  lead  to  a  fractional 
quantum  Hall  effect  because  interaction  energy  will  always  dom¬ 
inate  kinetic  energy  [142[.  To  date,  only  a  few  lattice  models  have 
been  proposed  which  are  expected  to  lead  to  a  fractional  quantum 
Hall  effect  [85,138-141  ].  The  relevant  figure  of  merit  in  such  models 
is  the  ratio  of  the  band  gap  to  the  bandwidth  of  the  fiat  band  with  a 
finite  Chern  number.  In  the  model  we  discuss  below,  we  find  this 
ratio  can  be  as  high  as  70,  which  is  among  the  largest  in  the  models 
reported  in  the  literature  thus  far  [85,138-141], 


Fig.  11.  Schematic  of  the  ruby  lattice  and  illustration  of  the  nearest-neighbor 
hopping,  t.t]  (real)  and  (complex),  and  the  three  types  of  “second-neighbor" 
spin-orbit  coupling  or  hopping  indicated  by  the  dashed  or  dot  dashed  lines,  t2.t3  and 
t4r  (from  Ref.  [85]).  (a)  The  spin-orbit  coupling  strength  t2  within  a  hexagon,  (b)  The 
spin-orbit  coupling  strength  t3  within  a  pentagon  composed  of  one  triangle  and  one 
square.  The  hoppings  t2  and  t3  are  present  on  all  such  bonds  of  the  type  shown  that 
are  consistent  with  the  symmetry  of  the  lattice,  (c)  Schemabc  of  the  hopping 
parameters  used  to  obtain  a  flat  band  with  a  finite  Chern  number  and  W/Eg  ^  70. 
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We  study  [85]  a  simple  s-band  tight-binding  model  on  the 
ruby  lattice  shown  in  Fig.  11.  For  fermions  with  spin,  the 
Hamiltonian  is  given  by  'H  =  Wruby+Wruby-so.  where 

'^ruby  =  -t  ^ 

ijeA.ff  A^A.u 

and 

Wruby-S0  =  !t2  ^  Vijf^pclcjp +  ((3  (19) 

on  the  ruby  lattice  shown  in  Fig.  11.  Here  cj^/Qa  is  the  creation/ 
annihilation  operator  of  an  electron  on  site  i  with  spin  cr.  As 
indicated  in  Fig.  ll(a,b),  t  and  H  are  real  first-neighbor  hopping 
parameters,  and  (2,  (3  are  real  second-neighbor  hoppings  (these 
appear  with  the  imaginary  number  i  in  Eq.  (19)  making  the  total 
second-neighbor  hopping  purely  imaginary  and  time-reversal 
symmetric).  The  quantity  Vy  is  equal  to  1  if  the  electron  makes 
a  left  turn  on  the  lattice  links  during  the  second-neighbor 
hopping,  and  is  equal  to  - 1  if  the  electron  make  a  right  turn 
during  that  process.  As  is  clear  from  Fig.  1 1 ,  the  unit  cell  of  the 
lattice  contains  six  sites,  so  six  two-fold  degenerate  bands  will 
result.  In  addition  to  the  real  hopping  parameters  t,t]  in  Eq.  (18), 
symmetry  also  allows  complex,  spin-dependent  nearest  neighbor 
hopping  with  imaginary  components  t',t]',  as  shown  in  Fig.  11. 
These  parameters  lead  to  a  variety  of  time-reversal  invariant 
phases  (including  Tl)  and  a  rather  complex  phase  diagram  [85]. 

Here  we  are  mainly  interested  in  the  spinless  case  with  broken 
time-reversal  symmetry  and  the  hopping  parameters  described  in 
Fig.  11(c).  For  spinless  fermions  we  study  Eq.  (18)  with 

t'  =  t+iffztj. 


t:'  =  tn  +  iffztu,  (20) 

and  add  the  hopping  inside  the  square  in  the  diagonal  directions, 
labeled  t4r.  This  is  shown  schematically  in  Eig.  11(c).  A  nonzero 
Chern  number  and  nearly  flat  band  occurs  in  this  case.  For 


r  M  K  r 


Fig.  12.  (a)  The  energy  bands  without  any  spin-orbit  coupling  and  with  it  tuned 
on,  given  by  Eq.  (18)  with  t  —  f  real  along  high  symmetry  directions.  Note  that 
there  is  a  Dirac  point  at  K  for  1/6  and  2/3  filling,  and  a  quadratic  band  touching 
point  at  r  for  1/2  and  5/6  filling.  (The  underlying  lattice  is  triangular,  as  it  is  for  the 
honeycomb  lattice.)  Note  also  the  flat  bands  along  the  T-M  direction  at  1/2  and 
5/6  filling,  (b)  The  energy  bands  with  t,  =  1.2t,tir  =  -l.Zt.ty  =  2.6t,t4r  =  — 1.2t  in 
Eq.  (20).  The  lowest  energy  band  has  Chern  number  —1.  The  ratio  of  the  band¬ 
width  to  band  gap  is  approximately  70  (from  Ref.  [85]). 


example,  with  t,  =  1.2t,fir  =  -1.2t,t]i  =  2.6t,t4r  =  -1.2t  the  Chern 
number  is  -1  and  the  gap  =  2.398t,  and  band  width=0.037t, 
which  gives  W/Eg^  70,  and  therefore  makes  an  excellent  candi¬ 
date  to  realize  a  fractional  topological  insulator  [85[.  The  corre¬ 
sponding  band  structure  is  shown  in  Fig.  12(b).  The  spinless  band 
structure  for  Eq.  (18)  with  t'  =  t  =  1  is  shown  in  Eig.  12(a).  Note 
the  interesting  band  crossing/touching  points  and  flat  portions  of 
the  bands  even  in  this  case. 

A  new  direction  has  emerged  in  the  search  for  systems 
with  flat  bands  with  conditions  favorable  to  realizing  fractional 
topological  insulators  or  fractional  quantum  anomalous  Hall 
states — transition  metal  oxide  interfaces  [33].  Research  in  this 
direction  is  rapidly  evolving  with  a  variety  of  interesting  interac¬ 
tion-driven  phases  possible  [30-32]. 

As  rich  as  the  physics  of  topological  insulators  and  quantum 
spin  liquids  is  in  two  dimensions,  the  behavior  is  even  richer  in 
three  dimensions.  This  is  because  the  three-dimensional  time- 
reversal  invariant  topological  insulators  are  qualitatively  different 
from  their  two-dimensional  counter  parts,  the  latter  of  which  can 
always  be  viewed  as  being  derived  from  an  underlying  quantum 
Hall  effect.  The  so-called  “strong  topological”  insulators  [143] 
(non-interacting)  are  an  entirely  distinct  three-dimensional  phase 
with  novel  response  properties  [8,51,143-148].  Moreover,  once 
interactions  are  brought  into  the  picture  the  physics  is  richer  still. 
We  now  turn  to  three  dimensions. 


4.  Three-dimensional  systems 

One  of  the  salient  features  of  Tls  is  that  their  boundaries 
possess  a  topologically  protected  “metallic”  state  that  is  robust  to 
disorder  [45,68,69,143,149-151].  Both  the  two  dimensional 
[116-123]  and  the  three  dimensional  [8,144-146]  boundaries 
have  been  shown  to  exhibit  interesting  responses  to  perturba¬ 
tions.  In  the  experimental  literature  on  Tls,  the  three  dimensional 
systems  have  dominated.  The  weakly  interacting  nature  of  topo¬ 
logical  insulators  has  enabled  accurate  predictions  based  on 
density  functional  theory  [78,147,152,153]  for  a  wide  range  of 
two  and  three  dimensional  systems  [138,154-167],  and  experi¬ 
ment  has  followed  with  confirming  data  in  a  large  and  rapidly 
growing  number  of  instances  [79,80,168-177].  There  are  also 
theoretical  predictions  for  systems  closely  related  to  Tls  with 
Weyl  fermions  in  their  bulk  and  unusual  boundary  excitations 
[32,178-182],  as  well  as  crystalline  topological  insulators  [183]. 

In  this  section  we  will  describe  systems  that  are  too  strongly 
interacting  to  allow  a  naive  application  of  density  functional 
theory.  On  the  most  exotic  side,  three-dimensional  versions  of 
exotic  fractional  topological  insulators  with  a  non-trivial  ground 
state  degeneracy  have  been  proposed  [184-187].  for  interacting 
systems  treated  at  the  mean-field  level,  a  spontaneous  generation 
of  spin-orbit  coupling  (with  topological  defects  possible)  can 
occur,  similar  to  two  dimensions  [125]. 

Here,  we  will  focus  on  semi-realistic  models  that  are  moti¬ 
vated  by  the  physics  of  transition  metal  oxides,  famous  in  recent 
years  for  high-temperature  superconductivity  [188,189]  and 
colossal  magnetoresistance  [190].  These  materials  typically 
involve  3d  orbitals  on  the  transition  metals  and  lead  to  strong 
Coulomb  interactions  due  to  their  localized  nature  (compared  to 
4d  and  5d  orbitals).  However,  4d  and  5d  orbitals  in  layered 
perovskites  such  as  Sr2Ru04,  Sr2Rh04,  Sr2lr04,  Na2lr03,  and  the 
hyperkagome  Na4lr308  are  more  spatially  extended  and  thus  the 
Coulomb  interaction  is  typically  weaker  than  those  with  3d 
orbitals  [191].  The  hyperkagome  Na4lr30g  is  an  important 
three-dimensional  quantum  spin  liquid  candidate  [192-196]. 

The  more  extended  nature  of  the  4d  and  5d  orbitals  compared 
to  the  3d  orbitals  leads  to  a  greater  level  splitting  in  a  crystal  field. 


CA.  Fiete  et  al.  /  Physica  E  44  (2012)  845-859 


855 


Fig.  13.  From  Ref.  [27],  (a)  An  illustration  of  the  pyrochlore  lattice  which  is 
composed  of  corner  sharing  tetrahedra.  Transition  elements  are  indicated  by  black 
solid  circles,  (b)  Each  transition  ion  is  surrounded  by  an  oxygen  octahedron  shown 
by  six  solid  blue  (dark  grey)  circles.  A  transition  ion  is  located  at  the  origin  of  the 
local  coordinate  and  is  shown  in  black.  We  study  a  trigonal  distortion  preserving 
C3  symmetry  applied  along  the  [111]  direction  (or  its  equivalent),  shown  by  two 
yellow  (grey)  faces,  and  an  elongation  preserving  C4  symmetry  along  the  z-axis  of 
the  local  coordinate,  (c)  A  schematic  representation  of  the  splitting  of  the  bare 
atomic  d-levels  (1),  due  to  a  cubic  crystal  field  arising  from  the  octahedral 
environment  (2),  unquenched  spin-orbit  coupling  in  the  t2g  manifold  (3),  and 
a  distortion  of  the  octahedron  (4).  The  values  of  the  splittings  in  (4)  depend  on  2 
and  ZI34. 

and  enhances  their  sensitivity  to  lattice  distortions  [197],  In  many 
oxides,  the  transition  ions  are  surrounded  by  an  octahedron  of 
oxygen  atoms,  MOe,  where  M  represents  a  transition  metal  ion. 
The  crystal  field  splits  the  five  degenerate  (neglecting  spin  for  the 
moment)  d-orbitals  into  two  manifolds  (see  Fig.  13c):  a  lower 
lying  tjg  (dxy.dyz.dzx)  manifold  and  a  higher  lying  eg  (d3z2_r2,dx2_y2 ) 
manifold  [198,199].  The  energy  separation  between  the  t2g  and  eg 
levels  is  conventionally  denoted  “lODq”  and  is  typically  on  the 
order  of  ~  l-4eV,  which  is  large  compared  to  many  3d  com¬ 
pounds  [200]. 

Besides  the  crystal  field,  the  relativistic  spin-orbit  coupling  is 
another  energy  scale  that  results  from  the  large  atomic  numbers 
of  heavy  transition  elements.  While  in  the  absence  of  spin-orbit 
coupling  the  on-site  Coulomb  interaction  is  of  the  same  order  as 
the  band  width  [201,202],  inclusion  of  strong  spin-orbit  coupling 
modifies  the  relative  energy  scales  [203].  Thus,  for  materials  with 
4d  and  particularly  5d  electrons,  one  expects  the  appearance  of 
novel  phases  with  unconventional  electronic  structure  and  order 
due  to  the  characteristic  energy  of  spin-orbit  coupling  approach¬ 
ing  that  of  the  Coulomb  interactions  [27,203,204]. 

In  this  section  we  focus  on  the  interplay  and  competition 
between  strong  correlation  effects,  spin-orbit  coupling,  and 
lattice  distortion  that  is  expected  to  be  important  in  heavy 
transition  metal  oxides  with  a  pyrochlore  lattice  [27].  In  the 
heavy  transition  metal  oxides  one  expects  both  the  spin-orbit 
coupling  [205,206]  and  the  lattice  distortion  energies  [197,207]  to 
be  of  the  order  of  0.05-0.5  eV,  while  the  interaction  energy  is 
typically  at  the  higher  end  of  this  scale  to  somewhat  larger,  0.5- 
2  eV  [178,179,205].  While  the  phase  diagram  of  an  interacting 
undistorted  pyrochlore  model  with  j  =  1  /2  has  already  been 
studied  [203],  we  expand  those  results  to  include  the  effects  of 
distortions  of  the  local  octahedra  on  the  phase  diagram  [27]. 

We  also  investigate  pyrochlore  oxides  at  different  d-level 
fillings  with  the  Fermi  energy  lying  in  the  quadruplet  j  =  3/2 
manifold,  which  has  not  been  considered  in  previous  works  [27]. 


One  of  our  motivations  is  to  see  if  the  j  =  3/2  manifold  can  also 
realize  the  interesting  Mott  phases  of  the  j  =  1/2  manifold  [203]. 
We  find  that,  indeed,  these  exotic  phases  can  be  realized  for  the 
j  =  3/2  manifold.  Moreover,  we  find  that  for  the  j  =  1/2  manifold 
“weak”  topological  variants  of  the  exotic  Mott  phases  can  also 
appear  in  the  phase  diagram  when  certain  types  of  lattice 
distortion  are  present. 

To  study  the  effects  of  lattice  deformations  [197,208],  we 
assume  that  the  octahedron  surrounding  an  ion  can  be  distorted 
in  two  ways:  (1)  a  trigonal  distortion  preserving  local  C3  sym¬ 
metry  and  (2)  an  elongation  (expansion)  of  octahedra  preserving 
local  C4  symmetry.  (See  Fig.  13b.)  The  former  has  been  argued  to 
be  rather  common  and  can  be  described  by  the  following 
Hamiltonian  on  each  transition  metal  ion  site  [197]: 

Tffri  =  —  ^^^{dy^dzx  -F  dy^dxy  -P  d^dxy)  -P  H.C.,  (21) 

where  4  3  parameterizes  the  strength  and  sign  of  the  C3  preser¬ 
ving  distortion,  and  the  C4  elongation/contraction  splitting  is 
described  by  [204] 

'Hel  =  A4IZ  =  A4(nyz-FTlzx),  (22) 

where  A4  parametrizes  the  strength  and  sign  of  the  distortion, 
and  Iz  is  the  z  component  of  the  effective  angular  momentum  of 
the  t2g  orbitals  related  to  the  occupation  of  the  d^y  orbital  by 
Uxy  =  ni—(lz)^  which  follows  from  the  constraint  =  rixy  +  nyz  +  rizx 
[204].  For  an  elongation  of  the  tetrahedron,  A4  <  0,  and  for  a 
compression  of  the  tetrahedron,  ZI4  >  0.  Trigonal  distortions 
appear  to  be  more  common  in  real  materials,  and  the  magnitude 
of  the  energy  splittings  can  be  crudely  estimated  from  density 
functional  theory  calculations  based  on  X-ray  determined  posi¬ 
tions  of  oxygen  atoms  around  the  transition  metals.  We  are  not 
aware  of  detailed  calculations  of  this  type  for  the  4d  and  5d 
pyrochlore  oxides,  but  closely  related  3d  systems  appear  to  have 
splittings  on  the  level  of  0.01-0.5  eV  [207].  We  take  this  as  crude 
estimate,  with  the  larger  end  of  the  energy  scale  probably  more 
likely  for  the  more  extended  4d  and  5d  orbitals. 

Thus,  the  local  Hamiltonian  describing  the  t2g  orbitals  on  each 
site  is 

'Hiocal  =  Hso+ntri  +  na,  (23) 

with 

Hso  =  -21  •  s,  (24) 

where  1=1  and  s=l/2  describe  the  orbital  and  spin  degrees  of 
freedom,  and  2  >  0  parameterizes  the  strength  of  the  spin-orbit 
coupling.  That  the  t2g  orbitals  can  be  effectively  described  by 
angular  momentum  1=1  comes  from  the  projection  of  the  d-orbital 
angular  momentum  into  the  local  basis  of  t2g  manifold  [203,204]. 

The  Hamiltonian  (23)  can  be  easily  diagonalized  and  its 
eigenvectors  describe  a  projection  onto  the  spin-orbit  plus 
distortion  basis.  We  will  denote  the  projection  by  a  matrix  M, 
which  contains  all  the  information  about  the  spin-orbit  coupling 
and  the  distortion  of  the  octahedra  (all  assumed  identical  so 
translational  invariance  is  preserved).  Moreover,  due  to  the 
presence  of  time-reversal  symmetry,  the  eigenvectors  form  Kramers 
pairs.  A  schematic  representation  of  splitting  t2g  upon  including  the 
terms  in  Eq.  (23)  is  shown  in  Fig.  13c. 

We  now  turn  to  a  derivation  of  the  effective  Hamiltonian.  We 
first  assume  A  =  A3  =  A4  =  0,  i.e.  neglect  the  contributions  in 
Eq.  (23).  To  obtain  the  kinetic  terms  of  the  Hamiltonian,  we  need 
to  describe  the  t2g  orbitals  of  a  single  ion  in  the  local  coordinate 
system  defined  by  the  octahedron  of  oxygen  atoms  surrounded 
the  ion,  and  we  need  the  p-orbitals  of  oxygen  in  the  global 
coordinate  system  [203].  The  hopping  of  electrons  from  one 
transition  metal  ion  to  a  nearest-neighbor  transition  metal  ion 
is  mediated  by  the  oxygen  p-orbitals.  (We  note  that  for  the 
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relatively  extended  5d  orbitals  direct  overlap  may  also  be  impor¬ 
tant,  as  well  as  further  neighbor  hopping  [209].)  We  thus 
compute  the  p-d  overlaps  to  determine  the  hopping  matrix 
elements.  The  local  and  global  axes  are  related  by  a  set  of  rotation 
matrices  [197,203].  The  combination  of  rotation  matrices  and  d-p 
overlaps  gives  rise  to  the  following  Hamiltonian: 

Hd  =  SaJ2dlyAy<y  +  t  Y.  (25) 

iy<7  ) 

where  /,  y,  and  cr  in  the  sums  run  over  lattice  sites,  t2g  orbitals 
(xy.yz.zx),  and  spin  degrees  of  freedom,  respectively.  The  8^  stands 
for  the  on-site  energy  of  the  degenerate  t2g  orbitals,  and 
f  =  ^pd7i/(Sp“Sii)  is  the  unrotated  hopping  amplitude  depending 
on  the  overlap  Integral  and  the  energy  difference  between  p 
and  d  orbitals.  The  parameter  t  sets  the  basic  hopping  energy 
scale  In  the  problem.  Without  loss  of  generality  we  set  =  0. 

The  effect  of  spin-orbit  coupling  and  distortion  are  included 
via  the  projection  of  Hamiltonian  in  Eq.  (25)  into  the  eigenvectors 
of  the  local  Hamiltonian  in  Eq.  (23)  using  matrix  M  as  follows: 

'^pyro  =  ^TgCiyCioi  -t  t  ^  •  (2®i 

ioc  (ia.i'a'} 

where  ua  (a=l,...,6)  stands  for  the  six  eigenvalues  of  local 
Hamiltonian  (23),  and  the  matrix  F  describes  the  hopping 
between  sites  given  in  the  local  basis  via  F  =  M*TlA.  The  cJ^(Cix) 
is  the  creation  (annihilation)  operator  of  an  electron  at  site  i  and 
in  local  state  a.  Finally,  we  add  a  Coulomb  interaction  to  obtain 


where  U  is  the  on-site  Coulomb  interaction  and  is  the  number 
of  electrons  on  the  5d  orbital  of  the  transition  metal  Ion.  In  the 
remainder  of  this  section,  we  investigate  the  zero-temperature 
phase  diagram  of  the  full  Hamiltonian  (27),  which  includes  the 
spin-orbit  coupling  and  lattice  distortions  in  Eq.  (23).  A  full 
description  of  our  results  is  given  in  Ref.  [27]. 

We  apply  the  slave-rotor  mean-field  theory  developed  by  Florens 
and  Georges  [25,26]  to  treat  the  effect  of  weak  to  Intermediate 
strength  Coulomb  Interactions  in  the  regime  where  the  charge 
fluctuations  remain  important.  In  this  theory  each  electron  operator 
Is  represented  in  terms  of  a  collective  phase,  conjugate  to  charge, 
called  a  rotor  and  an  auxiliary  fermion  called  a  splnon  as 

Qa  =  (28) 

where  Cjc.  is  the  electron  destruction  operator  at  site  /  with  quantum 
number  a,  representing  the  states  in  Eq.  (27).  The  factor  e‘®'  acting 
on  the  charge  sector  is  a  rotor  lowering  operator  (with  6,  a  bosonic 
field),  and  Is  the  fermionic  spinon  operator.  The  product  of  the 
two  results  in  an  object  with  fermi  statistics,  needed  for  the  electron. 
Note  the  rotor  part  only  carries  the  charge  degree  of  freedom  while 
the  splnon  part  carries  the  remaining  degrees  of  freedom  a.  There¬ 
fore,  an  electron  has  natural  spin-charge  separation  if  a  is  spin  in 
this  representation.  Substituting  Eq.  (28)  into  Eq.  (27)  and  applying 
mean-field  theory  results  in  the  phase  diagrams  shown  in  Figs.  14 
and  15  [27].  The  two  to  three  dimensional  crossover  and  issues 
associated  with  gauge  fluctuations  beyond  mean-field  theory  have 
been  discussed  in  Ref.  [210].  The  method  correctly  recovers  the 
H  ->  0  limit  in  mean-field  theory. 

There  are  several  features  that  are  Important  to  emphasize. 
First,  we  have  identified  a  new  phase — the  weak  topological  Mott 
Insulator  (WTMl) — In  the  regions  shown  in  Fig.  14(b)  ]27].  This 
phase  is  expected  to  have  gapless  thermal  transport  without 
charge  transport  along  a  certain  class  of  topological  defects,  as 
described  in  Ref.  ]28].  Also,  we  expect  a  rather  rich  phenomen¬ 
ology  In  the  spin  sector  of  WTMls,  similar  to  that  recently 
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Fig.  14.  Phase  diagram  of  the  j  =  1  /2-band  model  corresponding  to  nj  =  5  with 
positive  ZI3  =  2t  (upper  panel)  and  negative  trigonal  distortion  ZI3  =  -2t  (lower 
panel).  The  dashed  line  separates  the  rotor  condensed  phases  (below)  from  the 
uncondensed  phases  (above).  We  set  t=l,  and  the  phases  labeled  are  as  follows: 
strong  topological  insulator  (STl),  weak  topological  insulator  (WTl),  gapless  Mott 
insulator  (GMI),  topological  Mott  insulator  (TMl),  weak  topological  Mott  insulator 
(WTMl)  and  metallic  phases.  From  Ref.  [27]. 

emphasized  in  weak  topological  Insulators  (WTl)  [211].  Second, 
it  is  clear  from  Fig.  14  that  the  sign  of  the  trigonal  distortion 
makes  a  big  difference  in  the  complexity  of  the  phase  diagram, 
with  the  negative  trigonal  distortion  producing  a  more  interesting 
phase  diagram.  Third,  comparing  Figs.  14  and  15  one  sees  that  the 
j  =  1  /2  and  j  =  3/2  manifolds  behave  differently  for  the  same 
distortion.  Fourth,  the  lower  panel  of  Fig.  15  shows  the  phase 
diagram  in  the  distortion-spin-orbit  coupling  plane.  The  distor¬ 
tion  axis  can  be  viewed  as  a  pressure  axis.  Thus,  the  application  of 
pressure  to  transitional  metal  oxides  can  drive  the  system 
through  a  complex  set  of  phases,  some  of  which  are  topological. 
This  is  potentially  a  useful  experimental  parameter  when  search¬ 
ing  for  some  of  the  exotic  phases  predicted  here. 

Our  slave  rotor  results  make  a  natural  connection  to  quantum 
spin  liquids  in  three  dimensions  described  by  fermionic  excita¬ 
tions:  All  the  phases  above  the  dashed  lines  in  Figs.) 4  and  15  are 
time-reversal  invariant  spin  liquids.  These  phases  have  fully  gapped 
charge  degrees  of  freedom  and  spin  degrees  of  freedom  described 
by  the  fermionic  spinons  in  Eq.  (28).  At  the  mean-field  level, 
the  spinons  have  an  associated  band  structure  that  can  either  be 
gapped  and  topological  (like  a  Tl)  or  gapless  (like  a  metal).  Thus, 
the  spin  sector  has  an  intrinsic  single-particle-like  nature 
imposed  on  it  by  the  formalism.  Relaxing  this  assumption  is  one 
of  the  ingredients  needed  to  obtain  the  most  exotic  fractionalized 
topological  Insulators  in  three  dimensions  ]  184-1 87].  Indeed,  a 
fruitful  approach  to  studying  and/or  discovering  exotic  fractiona¬ 
lized  topological  insulators  in  three  dimensions  may  be  to  deepen 
our  understanding  of  quantum  spin  liquids  in  three  dimensions. 
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Fig.  15.  Phase  diagram  of  the  j  =  3/2  band  model  with  n^—3,  including  the 
trigonal  distortion  of  the  octahedra.  The  labeling  of  the  phases  is  the  same  as  that 
used  in  Fig.  14.  In  the  upper  panel  A3  =  2t.  In  the  lower  panel  the  interaction  is 
fixed  at  U  —  4t  and  the  strength  A3  of  the  trigonal  distortion  is  varied,  illustrating 
possible  phases  that  may  arise  upon  the  application  of  pressure  to  a  real  system. 
All  energies  are  expressed  in  units  of  t.  In  both  phase  diagrams  the  dashed  line 
separates  the  rotor  uncondensed  phase  (above)  from  the  condensed  phase 
(below).  We  note  the  “pocket”  of  STI  around  Arj  1,^3 1.5  in  the  lower  figure 
has  a  numerically  difficult  to  determine  boundary  with  the  metallic  phase;  we 
have  present  our  best  assessment.  From  Ref.  [27]. 

particularly  those  whose  spinous  are  strongly  correlated  and  not 
well  described  within  a  single-particle  formalism  [186].  In  a 
gapless  system,  these  would  be  “algebraic  spin  liquids”  [4],  One 
could  then  “add  back”  the  charge  degrees  of  freedom  to  obtain  an 
exotic  topological  insulator  [27],  However,  it  should  be  empha¬ 
sized  that  correspondences  between  TIs  and  QSLs  described  in 
this  paper  have  all  been  between  (JSLs  with  excitations  described 
by  fermions.  At  the  moment,  it  remains  unclear  how  to  relate  spin 
liquids  described  by  bosonic  excitations  (Schwinger  bosons)  to 
the  topological  insulators,  although  some  interesting  parallels 
have  been  noted  [212]. 


5.  Conclusions 

In  this  work,  we  have  emphasized  connections  between  topolo¬ 
gical  insulators  and  quantum  spin  liquids  from  a  variety  of  view¬ 
points,  with  an  emphasis  on  their  common  topological  structure. 
We  have  described  results  from  one,  two,  and  three  dimensional 
systems.  In  one-dimension  we  have  emphasized  the  important  role 
that  quantum  entanglement,  particularly  the  entanglement  spec¬ 
trum,  can  play  in  studying  and  classifying  phases. 

In  two  dimensions,  we  have  demonstrated  the  Kitaev- 
Haldane-Kane-Mele  correspondence  and  shown  how  it  can  play 
a  useful  role  in  finding  topological  insulators  and  interesting 
quantum  spin  liquids.  In  interacting  systems  we  have  described 


the  different  roles  Dirac  points  and  quadratic  band  crossing  points 
play  in  the  instability  to  spontaneously  generated  topological 
phases,  even  when  there  is  no  intrinsic  spin-orbit  coupling.  We 
have  given  an  example  of  a  lattice  model  of  spinless  fermions  that 
possesses  a  very  flat  band  with  a  finite  Chern  number.  This  model 
and  related  ones  can  be  used  to  help  guide  the  discovery  of 
fractional  topological  insulators  in  two  dimensions.  We  have  also 
emphasized  transition  metal  oxide  interfaces  as  a  promising 
experimental  system  for  the  realization  of  interaction-driven 
topological  phases. 

In  three  dimensions,  we  have  focused  on  semi-realistic  models 
appropriate  for  transition  metal  oxides  with  4d  and  5d  electrons. 
We  have  studied  the  interplay  of  distortion,  interactions,  and 
spin-orbit  coupling  and  presented  rich  phase  diagrams  for  dif¬ 
ferent  d-shell  fillings.  We  have  also  found  a  new  phase,  the  “weak 
topological  Mott  insulator”,  which  is  expected  to  reveal  itself  via 
topologically  protected  gapless  modes  along  certain  classes  of 
defects  [28[. 

Looking  ahead,  it  is  clear  that  two  of  the  chief  challenges  with 
three-dimensional  TIs  (which  are  most  likely  to  appear  in  appli¬ 
cations)  are  sample  quality  and  reliable  tunability  [29].  To  date, 
the  most  pressing  issue  is  to  obtain  samples  that  have  a  highly 
insulating  bulk  (so  that  most  electrical  conductance  is  along  the 
boundary).  In  recent  months  important  advances  have  been  made 
on  this  front.  Both  Bi2Te2Se  [213,214]  and  Bii.sSbo.sTei.ySei.a 
[215,216]  are  reported  to  have  more  than  70%  of  their  conductiv¬ 
ity  on  the  surface,  and  there  are  indications  that  TI  “work  horses” 
like  BbSea  are  improving  in  quality,  too  [217]. 

In  the  regime  of  very  strong  electron  interactions  (beyond 
where  the  slave-rotor  formalism  should  apply),  magnetic  phases 
with  broken  time-reversal  symmetry  become  more  likely  [203]. 
When  spin-orbit  coupling  is  strong  in  this  regime,  exotic,  con¬ 
ventionally  ordered  phases  are  also  expected  to  appear 
[195,196,204,218],  including  spin-orbital  liquid  states  [219,220]. 
It  seems  likely  that  the  new-found  capability  to  engineer  gauge 
fields  in  cold-atom  systems  will  play  a  role  in  our  understanding 
of  TIs  and  QSLs,  too  [221-226].  Many  novel  behaviors  certainly 
await  identification,  and  keeping  in  mind  the  “connections” 
between  various  topological  phases  seems  a  good  strategy  to 
help  guide  discovery. 
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